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I. 



I NTRODUCTigN 



For the Army to fight effectively in a resource scarce 
environment, the quantitative decision making techniques of 
operations research are important skills for Army staff 
officers. Staff officers are expected to be able to put 
numbers in their estimates when briefing commanders. They 
are expected to be able to measure and evaluate complex 
operations and subordinate units. They are expected to be 
frugal managers of time and money. And above all, staff 
officers must be able to apply sound, quantified reasoning 
in planning how to win the air-land battle. 

The use of hand-held programmable calculators by Army 
staff officers has the potential for improving the use of 
quantitative decision making techniques throughout the Army 
Faster and more accurate than paper and pencil, the calcula 
tor is less expensive and more portable than larger comput- 
ers. Sven when compared to the latest micro-computer 
systems or to portable terminals used for distributed data 
processing, the hand-held programmable calculator offers 
advantages in cost, reliability, power consumption and emis 
sion of electromagnetic radiation. Hand-held programmable 
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calculators have already been successfully used by soldiers 
in the field for applications in artillery fire direction, 
surveying, and navigation. In addition, large numbers of 
Army officers own their own pocket calculators and routinely 
use them for staff planning and reporting functions. 

In January of 1981 the 0. S. Army Command and General 
Staff College at Fort Leavenworth, Kansas selected a pro- 
grammable calculator for the Combined Arms and Services 
Staff School (CAS 3 .) Using both resident and non-resident 
instruction, this course is designed to teach all Army cap- 
tains staff techniques and procedures. As a significant 
part of the curriculum, the students are introduced to sub- 
jects such as statistics and regression, decision theory, 
combat modeling and linear programming. Considering the 
large number of officers projected to attend this course in 
future years, this course represents the most widespread 
training in operations research techniques ever attempted by 
the Army. The decision to provide a sophisticated calcula- 
tor to these students on an experimental basis was made for 
two fundamental reasons. First, the availability of a cal- 
culator with immediate field utility should motivate the 
student to apply the quantitative techniques as compared to 
the student who would be forced to do all calculations by 
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hand. second, the power of the calculator permits classroom 
discussion cf techniques such as linear programming and 
regression which are very difficult and time consuming to 
perform manually. 

This thesis documents the author's work to support the 
use of a calculator in the Combined Arms and Services Staff 
School. Initially, the intent was to produce a series of 
lesson materials incorporating the use of the calculator on 
a series of operations research topics which have immediate 
application for the Army division level staff officer. 
Instead, the work accomplished focused on the design and 
construction of a system to make the programming and testing 
of calculator programs easier and more efficient. Except 
for the introduction, this thesis is written for the person 
wishing to implement the programming support system 
described. The implementor must have a detailed knowledge 
of the instruction set and programming characteristics of 
the HP41CV calculator as described in Wickes [Ref. 1: pp. 

6-20]. For the eventual user of the system, as compared *o 
the implementor, the system itself provides on-line documen- 
tation on how to use the system and what commands and 
options are available. Figure 1 shows the command menu dis- 
played on the terminal screen by this interactive program; 
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Figure 2 gives a more detailed explanation of each of the 
commands; and Figure 3 displays the on line introductory 
material that is provided to new users of the system. For 
the user, a knowledge of the information contained in the 
calculator owner's handbook [Ref. 2 ] is sufficient to begin 
writing calculator programs using the support system 
described. 

The calculator selected by the Command and General Staff 
College, the Hewlett-Packard HP41CV, typifies the state of 
the art in off-the-shelf calculator technology. While not 
without disadvantages, this calculator was selected because 
of its power and features which make it easier for Army 
staff officers to use. First and most important of these 
features is the ability of the calculator to manipulate 
alphabetic characters in addition to numeric data. The cal- 
culator can display the name of a variable when input data 
is required or label output when the calculation is com- 
pleted. With this feature, the calculator helps the user 
know what data to input or what action to take next. It 
also helps alleviate the need for constant reference to 
printed instructions which are difficult to use under field 
conditions. 
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Figure 1: Programming Environment Command Menu 
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Figure 2: List of Commands 
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Figure 2 (Continued) 
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Figure 3: On-Line Introductory Material 



A second important feature is the multiplicity of means 
by which programs can be entered inro the calculator. Mag- 
netic cards, read only memory, and optical bar code are all 
available and each has advantages depending on the situa- 
tion. For the long term, read only memory offers the abil- 
ity to retain very large programs (in excess of 8000 bytes) 
and the simplest and most reliable means of entering pro- 
grams into the calculator under field conditions. For the 
short term, optical bar code offers the least expensive 
method of reproducing and distributing calculator software 
that has not been subject to extensive field testing. In 
addition, as shown in this thesis, the optical bar code can 
provide an important link between a main-frame computer and 
the hand-held calculator. 

A third important feature of the HP41CV is its rela- 
tively large memory capacity as compared to programmable 
calculators such as the Texas Instruments TI-59. A large 
amount of memory permits the solution of larger, often more 
realistic problems than could previously be solved on a 
hand-held device. A demonstration program given in this 
thesis for linear programming is an example of an applica- 
tion where the full memory capability of the HP41CV is 
required :o be able tc solve realistic problems. 
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To take advantage of the calculator's unequalled economy 
and portability, the operations research analyst is chal- 
lenged to overcome its limits of speed and memory capacity. 
The preparation of calculator software is as difficult, if 
not more so, than the preparation of software for larger 
computers. To accomplish the most possible with the hand- 
held device, the calculator programmer is often forced to 
write programs which are very difficult to comprehend when 
examined by other programmers. As Dahl, Dijkstra and Hoare 
[Ref. 3: pp. 1-10] point out there are limits to human com- 

petence which interfere with the programming process. In 
the past, with less mature calculators which constrained the 
typical program to a few hundred program steps, these limits 
to human competence were neither as apparent nor as economi- 
cally important as they are with the HP41CV. Accordingly, 
it is not envisioned that the average Army officer who uses 
the HP41C7 on real world problems which push the calculator 
to the limits of its capability would write their own pro- 
grams. In particular, it was never intended that the stu- 
dents in the Combined Arms and Services Staff School would 
be taught calculator programming. It is a tribute to the 
power of the device and the quality of the calculator soft- 
ware when a relatively inexperienced user can run complex 
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programs using little more than the digit entry keys and the 
run-stop key on the calculator. This does not mean that the 
user must not have a clear understanding of his problem or 
the solution technique, but rather it means that the calcu- 
lator should not require programming skill or extensive 
training prior to application. 

The growing complexity of calculator programs described 
above and the realization that calculator programs for Army 
field use are not programmed in the field, suggest the need 
for a system to support the development, distribution and 
maintenance of calculator software. An operations research 
analyst or other professional programmer must be able to 
more efficiently prepare calculator programs than by keying 
them into the hand-held device. By preparing the programs 
initially on a larger computer, such as the IBM 3033, the 
programmer can use the speed and storage capability of the 
larger machine to great advantage. In addition, the avail- 
ability cf a full-screen video text editor speeds the pro- 
cess of program revision and maintenance. By providing a 
capability to integrate comments directly into the source 
code on the larger computer, program documentation is more 
easily provided. Essentially the idea is that a programmer 
would write the calculator program using a terminal 
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connected to a large computer. After the calculator program 
is entered into the large computer, a compiler program 
running on the large computer would check the calculator 
program for errors and convert the mnemonic instructions 
into the "key codes" which are the numeric instructions 
actually executed by the calculator. Then an em ulator 
program running on the large computer would take the numeric 
instructions from the compiler and execute the program--in 
effect making the large computer produce the same effects as 
the calculator only much faster and more efficiently for the 
programmer. Finally, when the program has been written and 
tested on the large computer, optical bar code is produced 
which allows for the economical distribution and use of the 
program in the field. To encourage the calculator program- 
mer to use the system described, this process should occur 
in an interactive programming environment in which the user 
can move from one step to another by issuing simple commands 
such as those listed and described in Figure 2 and receive 
help or cn line instruction whenever desired. Under this 
proposed system, the advantages of both the larger computer 
and the hand-held calculator are used appropriately in a 
mutually supporting manner. This thesis presents two of rhe 
components of this proposed sys-em. First, an I3M EXEC II 
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program is given which provides an interactive programming 
environment for users operating under IBM's Conversational 
Monitor System (CMS.) A short discussion of the design of 
this program and a complete copy of the source code is 
contained in Appendix D to this thesis. Secondly, a cross 
compiler written in IBM standard FORTRAN IV is provided for 
translating calculator mnemonic instructions into the key 
codes necessary for use by the emulator and also for the 
production of optical bar code. The term cross compiler 
refers to the fact that the program runs on one machine (the 
larger computer) but compiles programs for another machine 
(the calculator.) A discussion of the design of this 
program and a complete copy of the source code is contained 
in Appendix D. To make the program easier to understand and 
adapt to new requirements, it is modularized into 24 
subroutines and is heavily commented. 

To illustrate the use of the system, two of the six 
example programs originally planned are provided in this 
thesis. Revised plans now call for the remaining four exam- 
ple programs to be issued at a later date as Naval Postgrad- 
uate School technical reports. Because the reasons for the 
delay constitute some of the most important lessons learned 
from this thesis research. Chapter 2 documents the process 
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with a technical discussion of the factors involved. The 
major conclusions described in Chapter 2 are the need for a 
prioritized list of criteria with which to evaluate calcula- 
tor programs and the need for more structure in the program- 
ming process. Chapter 2 is technically oriented and assumes 
the reader is familiar with the concepts of structured 
programming . 

Each of the calculator program examples is described in 
a separate appendix in which the documentation listed in 
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Figure 4: Components of Program Documentation 

Figure 4 is provided. The first example on single variable 
statistics is documented in Appendix A and uses the calcula 
tor in an area where calculators have long been used, but 
does so in a way that shows the unique capabilities of the 
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HP41CV. A second example on linear programming is docu- 
mented in Appendix B and illustrates an area where calcula- 
tors have not received widespread application. Most 
calculator linear programs which have been published to date 
have been either incomplete algorithms or have been limited 
to very simple problems. 

A third example, which by its nature does not conform to 
the documentation standards outlined above, describes a set 
of utility routines which could be distributed in read only 
memory. Programs for read only memory have different char- 
acteristics from other calculator programs and Appendix C is 
provided to illustrate some of these differences. 
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II. THE PROGRAMMING ENVIRONMENT 
A. CHAPTER OVERVIEW 

This Chapter examines calculator programming within rhe 
context of the author's experience in preparing HP41CV pro- 
grams in support of the Combined Arms and Services Staff 
School. With the advanced capabilities and features of the 
HP41CV# it was hoped that a complete package of software 
could be prepared quickly. To document why this did not 
occur# this chapter will examine strengths and weaknesses of 
the calculator in relationship to a collection of techniques 
referred to in computer science as structured programming. 
For the reader unfamiliar with this term# the previously 
cited work by Dahl# Dijkstra, and Hoare [Ref. 3] is 
recommended. This chapter is technically oriented and does 
assume familiaritiy with structured programming concepts. 

When programming calculator programs for personal use, 
most programmers, including the author, do nor find the task 
difficult. Programming a hand-held calculator with the 
capabilities and features of the HP41CV can be a rewarding 
experience. It is rewarding to master rhe algorithm of an 
operations research technique on a hand-held device. The 
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educational value in programming the calculator has been 
recognized by many educators, including Hamming [Ref. 4: 
pp. 2-3] and Weir [Ref. 5: pp. xii-xiii]. Providing a pro- 

gram for general distribution which makes optimum use of the 
calculator is quite a different situation. It was the 
author's experience that programs, which gave correct 
answers when used by the author, often had to be completely 
re-written several times before being acceptable. This 
problem became more acute as the size of the programs grew 
beyond 400 program steps, for at that size it became 
increasingly difficult to modify programs without affecting 
the total design. The major conclusions described in this 
chapter are the need for a prioritized list of criteria with 
which to evaluate calculator programs and the need for more 
structure in the programming process. 

B. STRUCTURED PROGRAMMING WITH THE HP41CV 
1 . The Need fo r S tru ct ur e 

To increase the efficiency of the programming pro- 
cess, a collection of techniques known as structured pro- 
gramming has received widespread attention in the computer 
science community. While there is no one definition cf 
structured programming, it does require three essential 
characteristics. First, there must be a logical structure 
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to the program which reflects the nature of the problem to 
be solved and any constraints imposed upon the solution. 
Second, the systematic process of stepwise refinement is 
used to limit the complexity of program segments. Third, 
the programming language must reflect the logical structure 
of the program and assist in stepwise refinement. These 
three characteristics represent not so much a detailed rec- 
ipe for program development as they do a philosophy of how 
programs can be more efficiently written. It was with this 
philosophy in mind, that a calculator programming support 
system was proposed which could take into account the 
strengths and weaknesses of the calculator; balance the 
structured programming philosophy with the other criteria 
listed below; and thereby solve the problems encountered in 
writing calculator software for the Combined Arms and 
Services Staff School. 

2. F un dame ntal Limitations of Calculato rs 

Writing programs to solve complex problems on a 
hand-held calculator is difficult both because of inherent 
limitations in the calculation speed and memory capacity of 
the machine and also the inability of the calculator's 
native programming language to directly support structured 
programming constructs. In many respects, the task is 
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similar to writing assembly level language programs for 
larger computers. Calculator programming features a 
powerful instruction set including advanced mathematical 
functions but lacks any ability to refer to variables by 
name instead of storage address. Like assembly language, 
the calculator's programming language consists of short 
mnemonic instructions typically followed by the storage 
location of the data to which the operation is to be 
applied. While a large amount of computer programming is 
still done in assembly language, it is generally accepted 
that programming in a higher level language such as FORTRAN 
is preferable. Programs written in an assembly language 
take more time to write and are not as easily changed as 
higher level language programs. Also, because they depend 
on the instruction set of a particular machine, they can not 
be easily transfered from one computer to another. These 
same disadvantages apply to calculator programming. In 
addition, because the hand-held device does not have the 
speed and memory capability of the larger machine, the 
calculator programmer must be even more mindful of the need 
to optimize his program to save program steps and execution 
time . 
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3. Modular Design 

The HP41CV supports structured programming as well 
or better than any other hand-held calculator. As described 
in the owner’s manual [Ref. 2: pp. 177-196], the machine 

primitive instruction XEQ encourages the construction of 
modular programs using calculator subroutines. Each subrou- 
tine can be a self-contained unit capable of being written 
and tested independently and used by multiple programs. 

This modularity is most strongly encouraged when routines in 
read only memory are used, for then the application program- 
mer can significantly reduce the number of program steps in 
his own program. This modularity, however, is not complete, 
since all variables are globally referenced and can be 
changed deliberately or inadvertently by any subroutine., 
This problem is no more apparent than with the use of read 
only memory, since one of the most limiting factors in using 
the read only memory programs as subroutines is conflict in 
the use of common registers. Also, unlike the modularity 
required in truly structured programs, there is no 
restriction limiting a subroutine to a single entry and a 
single exit point. In structured programs, such limits on 
entry and exit serve to define the fundamental building 
blocks by which stepwise refinement is made possible. With 
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the calculator, however, multiple entry and exit points are 
most useful for allowing a common routine to handle a 
duplicity of problem conditions. In this thesis, for 
example, programs are given for which two standard entry 
points are provided. One entry point uses an alpha-numeric 
label and an audio prompt to speed data entry, while a 
second entry point uses the alpha-numeric label but 
suppresses the audio tone. After data entry, the value 
entered is displayed, and the user is required to verify the 
accuracy of the data entered. By using the same subroutine 
with different entry points, memory space is saved overall 
at the sacrifice of the structured programming philosophy. 

4 . C ontrol of Pr o gra m Flow 

A basic deficiency prohibiting the HP41CV from 
directly supporting structured programming is the way in 
which program flow is controlled. Programming languages 
which support structured programming typically have instruc- 
tion constructs such as WHILE — ENDWHILE, REPEAT-- UNTIL, or 
LOOP — QUIT--ENDLOOP which make programming loops clear and 
concise. Constructs such as IF — THEN — 2LSEIF--ELSZ--ENDIF 
and the CASE statement make the evaluation of conditional 
expressions efficient and relatively error free. Also, 
structured programming languages typically discourage the 
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use of GOTO unconditional transfers because they lead to 
confusing code. In contrast, the HP41CV programmer must 
write his own looping constructs and his own conditional 
evaluation constructs using machine primitive instructions 
which somewhat obscure the program's basic objective and 
flow of control. In addition, it is difficult to avoid dis 
turbing pending operations in the stack registers when a 
conditional statement must be evaluated. As can be seen by 
the short program shown in Figure 1, the notation of the 
programming language does not permit structured program 
flow . 

5. C larifica tion of Program Structure 

Because no calculator, including the HP41CV, sup- 
ports named variables, the use of comments as an integral 
part of the calculator program is vital if the logical 
structure of the program is to be made clear as required by 
structured programming. Comments should provide the vari- 
able names when storing and recalling data; they should pro 
vide clarification of program flow; and they should mark 
subroutine boundaries and entry and exit points to make it 
easier to identify segments of the program. With the 
HP41CV' s stack oriented architecture, it is also frequently 
useful to display the names of the contents of each of the 
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Given rhe number n in the x-register, this program 
fragment will sum the data values stored in memory 
locations 1 through n. 



Instruction 


Comment 


LBL "SUM 
1E3 

/ 

1 


To execute press "XEQ SOW". 


+ 


Establishes a loop counter. 


0 


Clears x and pushes loop 


LBL 0 0 


counter into y. 


RCL IND Y 


Recall the next data value. 


+ 


Accumulate the sum. 


ISG I 


Increment the loop counter. 


GTO 0 0 


If more data remains, branch; 


RTN 


else, quit and display sum. 



Figure 5: Example Program to Add n Numbers 
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stack registers. In Appendix C on common subroutines with 
read only memory application, a shell sort [Ref. 6: pp. 

84-95] routine is given which employs the technique of using 
comments to display the names of the variables on the stack 
register . 

6. DalS T y pes an d Indirect Addressing 

Calculator programs represent more than a sequence 
of keystrokes; they also represent the manipulation and 
transformation of data. For maximum efficiency, the manipu- 
lation of data should be structured so as to prevent common 
programming errors. For this reason, most computer lan- 
guages which directly support structured programming enforce 
data type correspondence between data and operations. Fre- 
quently the formal declaration and initialization of vari- 
ables is also required. The HP41CV handles two types of 
data — real numbers and alphanumeric characters. While no 
formal declaration of variables is required, type checking 
is done automatically and is transparent to the user. Any 
attempt to perform an arithmetic operation on alpha-numeric 
data will result in the message "ALPHA DATA" and the program 
will halt. 

Because there is no formal declaration of variables, 
the programmer writing programs for the HP41CV must use 
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extreme caution in managing his dara set and insuring that 
the numbers stored and recalled by the calculator program 
are in fact the data elements desired. A typical example of 
an improper data reference occurs when a program is using 
indirect addressing and attempts to store or recall data 
from a non-existent data register. This programming error 
is so common that a special error message "NONEXISTENT” is 
provided by the calculator when this error is detected. 
Indirect addressing is an important feature which gives the 
calculator a considerable amount of power and flexibility, 
but also represents an additional responsibility for the 
programmer to explicitly control. On the HP4 lev all indi- 
rect addressing calculations must be specifically provided 
by the application progra m- -the re are no vector or array 
data types such as usually found with higher level lan- 
guages. In an attempt to make indirect addressing more 
transparent to the programmer, an experimental subroutine 
was prepared to recall an arbitrary element of a matrix 
stored as a two dimensional array. This subroutine, which 
is shown in Figure 2, was used in a simultaneous differen- 
tial equation combat model and the results evaluated. It 
accomplished the task, but slowed the execution of the pro- 
gram considerably (resulting in an overhead of 10.5 seconds 
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of extra execution time for every 100 subroutine calls) ana 
did not significantly improve the size or legibility of the 
application program. Accordingly, this technique is not 
recommended and indirect addressing remains a task that must 
be treated explicitly by the application programmer. 

C. ADDITIONAL CRITERIA FOR PROGRAM EVALUATION 

Calculator programming in many respects resembles a 
multi-criteria decision problem. On the surface the 
criteria for program effectiveness are quite straight 
forvard--the program must yield the correct answer, run 
quickly, require the fewest possible memory registers and be 
user friendly. Unfortunately, these objectives often 
conflict and can not always be simultaneously achieved. In 
particular, the principles of structured programming are 
often in conflict with the desire to reduce the size of 
programs and increase their execution speed. It is also 
true that the objectives of structured programming concern 
the process of writing programs, whereas the additional 
criteria listed concern the final program product itself and 
are therefore logically considered separately. Attempting 
to achieve all criteria at once can lead to failure, and 
some tradeoffs must be considered to evaluate programs and 
guide program development. The following criteria represent 
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Entry to this routine assumes the x register contains 
the column number and the y register contains the row 



number. The base address must be stored in R04 and 
the dimension of the matrix must be stored in R05. 






1 LBL "RCLM 

2 RCL 04 

3 + 

4 XOY 

5 1 

6 - 

7 RCL 05 

8 * 

9 + 

10 RCL IND X 

11 RTN 

12 END 



(BASE ADDRESS REGISTER 
(ADD BASE TO COLUMN NUMBER 
(RECALL THE ROW NUMBER 



(DIMENSION OF THE MATRIX 

(ADDRESS IS NOW IN X REG 
(RECALL THE DATA DESIRED 



Figure 6: Program to Recall an Element of a Matrix 
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"lessons learned” in developing application programs as 
examples for this thesis. 

1 • User Fri endline ss 

User friendly programs consider the application 
environment and do not task the user to be all knowing or 
without error in entering data. While individuals differ 
greatly with experience, the average user will make frequent 
errors in entering data with the hand-held calculator's 
small keyboard. In talking with officers who had used the 
TI-59 calculator in the field for fire direction, it was 
discovered that most preferred to use the printer with the 
calculator because it allowed data to be checked after 
entry. This was in spite of the fact that the printer and 
calculator combination is more costly, less portable and 
less suitable for use in the field than the calculator 
alone. In short, user friendliness was more important than 
these other criteria. For this reason, it should be manda- 
tory that any calculator programs intended for Army use in 
the field must allow the verification of data after entry. 
Because the use of the printer obviates many of the advan- 
tages of the hand-held calculator, the printer should not be 
required for this verification. Dne of the considerable 
advantages of the HP41CV is that the large amount of program 
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memory makes it possible to store the input values and 
perform this verification. However, programs written with 
this criteria in mind may not appear be most efficient to 
the casual observer. 

Another important aspect of user friendliness is 
limiting the complexity of the calculator and the actions 
reguired to get results. The typical Army officer has lit- 
tle appreciation for the multitude of scientific and mathe- 
matical functions labeling the keys of the HP41CV. Yet the 
common programming practice of using the toe two rows of 
calculator keys to indicate the identity of a variable 
either upon input or output increases confusion over the use 
of the function keys. This works as follows: When local 

alphabetic labels are used in a program to represent entry 
points by which a user indicates the identity of an input 
variable or reguests a particular output variable, then the 
first two rows of keys on the HP41CV become subroutine exe- 
cution keys pointing to these local labels when the calcula- 
tor is in user mode. This feature was very important on the 
HP67 and TI-59 where the lack of alpha-numeric capability 
reguired this method of program execution in order to most 
easily determine the identity of the input or output value, 
but it is less important on the HP41CV. It is almost always 
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true that a program which requires the use of local labels 
is harder to use, and requires more frequent reference to 
the user instructions than a program which uses only the 
run-stop key and properly prompts the user and labels output 
values. 

2. Ex ecution Speed 

The second most important criteria for a calculator- 
program is that it must yield results relatively quickly. 

In preparing example programs for this thesis, this point 
became very clear when testing two particular programs. One 
program, a simultaneous differential equation combat model, 
required in excess of 150 data values in order to yield 
results. It should be noted that it was only with the 
introduction of the HP41CV that it became feasible to con- 
sider such large problems on a hand-held device. To accomo- 
date the size of the model, the program was written so as to 
economize on program steps at the expense of increased exe- 
cution time. It became immediately obvious upon initial 
testing that this had been the wrong priority--f or users of 
the program were not impressed with either the use of the 
calculator or the utility of the combat model. If such user 
acceptance is not present, then the calculator program will 
remain unused, no matter how elegant the design to conserve 
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memory. In contrast, the linear programming example given 
in Appendix B was written so as to emphasize speed even if 
it meant including code redundancy. This program has been 
well received in part because it is so much faster than 
paper and pencil methods. 

The easiest and most effective technique that is 
useful in increasing speed is to decrease the number cf pro- 
gram steps that the calculator must process inside program 
locps. For example, if two different program options 
require similar but slightly different actions within a pro- 
gram loop, it is tempting to insert a program flag check and 
branching instructions within a loop so as to use the same 
loop for both conditions. But this means that the calcula- 
tor must test the flag and branch inside the loop even 
though the program is probably shorter overall. 1 

Instead, if the application permits, the memory capacity 
of the HP4 1CV can be used to best advantage by testing the 
flag once and then providing separate program loops for the 
two conditions. Again, this does not appear elegant to the 
casual observer, but it may result in a more successful 
program overall. This principle was discovered while 

1 Branching is required when the flag tests either set or 
clear if more than one instruction is required to account 
for the differences in the two conditions. 
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programming the single variable statistics program given in 
Appendix A. Initially, this program used a common loop for 
all data input and output operations, including reviewing 
the input data and making individual corrections. By 
providing a separate, somewhat redundant loop for data 
correction, the time required to input data points was 
reduced . 

D. A PROGRAMMING SUPPORT SYSTEM 

Considering the structured programming philosophy dis- 
cussed above in paragraph B and the additional criteria for 
evaluating programs listed in paragraph C, it becomes imme- 
diately obvious that programming with the calculator alone 
will never meet even a majority of these objectives. It 
must be recognized that the problem under consideration is 
not how the average person who owns a calculator should pro- 
ceed to program it for his own personal use, but rather how 
the Army can best provide the most cost-effective computa- 
tional resource for field use. For these reasons, a 
comprehensive programming support system is required. The 
programming support system outlined here will consider only 
the requirement for cost-effective preparation and 
maintenance of the calculator programs and not the broader 
issues of distribution and logistic support for the entire 
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calculator system to include hardware, training materials 
and printed references. 

1. A Cross C ompiler and 3 a r Code Generator 

The first requirement for an operational support 
package is to free the programmer from the limitations of 
the hand-held calculator itself. Sven with the printer and 
other peripherals, the calculator is no match for the larger 
machine when large programs must be examined or edited. In 
addition, the calculator is not currently capable of produc- 
ing its own optical bar code as required for economic repro- 
duction and distribution of the software. Accordingly, a 
cross compiler for the HP41CV was listed as the first 
requirement of the programming support system. Such a cross 
compiler has been written and is the major outcome of this 
thesis effort. This cross compiler accepts an HP41CV pro- 
gram written in the language of the calculator and returns 
the finished bar code as output. Any valid HP41CV program 
will be processed without need for modification by the cross 
compiler. In addition to the basic language of the calcula- 
tor, the user is allowed to inject comments directly into 
the source code with the use of the left parenthesis as a 
comment indicator mark. The ability to make comments 
directly in the source code makes the calculator programs 



40 



more legiole and more easily modified at a later date or by 
another programmer. Often, well placed comments can make up 
for a lack of structure in the program itself as far as 
legibility and maintainability are concerned. Having the 
comments directly in the source code facilitates their use 
and helps insure that they are as up to date as the program. 
For the average programmer, use of unmodified HP41CV source 
code augmented with a comment indicaror will represent the 
most common use of the cross compiler. The cross compiler 
is described in more detail in Appendix D including a 
complete listing of the source code. 

2. A Ca lcula tor Emulator 

After the calculator source code has been processed 
by the cross compiler, a need exists to be able to run the 
program without the wait for the generation of bar code. In 
addition, for the future development of read only memories 
for the calculator, an emulator program is required because 
the calculator itself can store only up to 2000 instructions 
in active random access memory. The read only memory can 
store up to four times this amount. Thus, the calculator by 
itself may not be capable of testing extremely large pro- 
grams or programs with large amounts of constant data also 
stored in the read only memory. Although an emulator was 
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not written for this thesis, the design of the cross 
compiler reflects the need for such a program. For example, 
the cross compiler generates an intermediate array of deci- 
mal integers which repesent the machine language of the 
HP 4 1 CV prior to conversion to binary. It was intended that 
these decimal integers could be used without modification or 
further translation within a FORTRAN computed goto state- 
ment. Thus, with the difficult translation, instruction 
parsing and syntax recognition already performed by the 
cross compiler routines, the emulator could consist of one 
large FORTRAN loop wherein a decimal integer was addressed 
in the instruction array by a program pointer variable. The 
integer is then immediately sent to a computed goto state- 
ment which would branch to the appropriate line of FORTRAN 
code which would simulate the referenced instruction, 
including updating the stack and the program pointer as 
appropriate. 

3. A Higher Level Lang uag e Compiler 

The final component in the calculator programming 
support system would be a program that would translate a 
higher level language such as PASCAL into HP41CV language 
which could then be sent to the cross compiler for verifica- 
tion and generation of the bar code and intermediate 
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calculator language listings. It is the higher level 
language compiler that would most directly make up for the 
weakness of the calculator in supporting structured program- 
ming. If would be able to increase the modularity of pro- 
grams, provide for named variables, make indirect addressing 
transparent and provide structured statements such as 
WHILE — ENDWHILE and IF— THEN — ELSE. Again, the design of 
the cross compiler anticipates this requirement and provides 
a considerable number of subroutines that would also be 
required by a higher level language compiler. These subrou- 
tines include a complete set of string functions for manipu- 
lating character data in FORTRAN and an instruction parser. 
3ecause it was envisioned that the higher level language 
compiler would also be able to process statements entered 
directly as HP41CV instructions, the cross compiler is con- 
structed so that the routine which compiles individual lines 
of HP41CV source code could be called as a subroutine by the 
higher level language compiler. Thus, all three major com- 
ponents of the proposed calculator programming support 
system would work together efficiently. 
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APPENDIX A 



SINGLE VARIABLE STATISTICS EXAMPLE 

INTRODUCTION: 

Calculating single variable statistics is one of the 
most frequently used applications of programmable calcula- 
tors. Army division level staff officers use single vari- 
able statistics to summarize and describe data for command 
briefings and periodic reports. The text by Mendenhall, 
Scheaffer and Wackerly [Ref. 7: pp. 3-13] is recommended as 
an introduction to the statistical measures calculated by 
the program given in this appendix. This program 
automatically calculates: 

• Mean and Median 

• Sample Standard Deviation 

• Sum of the Squared Deviations about the Mean 

• Coefficients of Skewness and Kurtosis 

• Minimum, Maximum and Range 

• Histogram Cell Frequencies 
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A single variable statistics program has been given as 
an example because of its immediate utility to the staff 
officer and to illustrate several features of the HP41CV 
which make it a superior device for Army field use. The 
most important of the these features is alphanumeric prompt- 
ing for input data values. The program given in this 
appendix provides an alphanumeric prompt for every input and 
output value and requires only the digit entry keys and 
run/stop key for data entry. Another important feature of 
the HP 4 1CV used by this program is its large memory 
capacity. This program retains up to 219 data points in the 
calculator's memory to allow the user to review the input 
data and make corrections during data entry. The large 
amount of memory allows the calculator to sort the data and 
calculate the order statistics including the minimum, 
maximum and median. Calculation of the median is a feature 
of this program which distinguishes it from other calculator 
statistics programs. In addition, without having to 
re-enter the data, the histogram may be calculated with a 
varying number of cells or a varying cell width. 

PROGRAM DESCRIPTION: 

The single variable statistics program has entry points 
for two different techniques of data input. The fastest 
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method, which provides both an alphanumeric prompt and an 
audio tone to speed data entry, may be called by execution 
of the program from entry point "STAT 1 . " A slower method, 
which provides greater accu racy and suppresses the audio 
tone for classroom use, may be called by execution of the 
program from entry point "S 1 . " When called from "SI," the 
program requires the verification of each data point after 



entry. The sequence of actions is as follows: 



1. The calculator displays an alphanumeric prompt. As 
an example, "XI?" is the prompt for the first point. 

2. The user enters the data value with the digit entry 
keys and presses the run/stop key. 

3. The calculator displays the data entered with a label 
derived from the alphanumeric prompt. For example, 

"x 1=3. 141 5" is a typical calculator response. This 
display is prompting the user to verify the correct- 
ness of the data displayed. 

4. If the value is correct, then the user simply presses 
the run/stop key and the calculator advances to the 
next point. 

5. If the value is erroneous, the user enters the correc 
value with the digit entry keys and then cresses the 
run/stop key. Then the calculator will again repeat 
srep 3 and ask the user to verify the data value. 

This process will continue until the user makes no 
modification to the data value. 

To run the program from either entry point the user may use 



the XEQ key, or assign the entry point label (" STAT 1 " or 



"SI") to a key and execute it by pressing that key in the 



USER mode. Further instructions on running programs and 
making key assignments are contained in the calculator 
owner’s manual [Ref. 2: pp 114-116], 
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In addition to the two initial entry points described 
above, several other alphabetic labels provide the user with 
functions that are called outside the normal sequence of 
program execution. Label " SR" provides the user with the 
capability to review the data stored in calculator memory, 
either before or after the data has been sorted. When used 
before the sort, the "SR" function is most useful in verify- 
ing the entire data set at one time. If used for this pur- 
pose, it should be called after all of the data has been 
entered and the mean of the data set is displayed with the 
"X BAR" label. If flag 21, the printer enable flag, is set 
"on" during this data review, then the calculator will stop 
as each point is displayed and the user may make corrections 
in the same manner as described above for the point-by-point 
verification associated with the "Si" entry point. When 
used after the sort, the "SR" function is most useful for 
displaying the order statistics for the data set. If used 
for this purpose, it should be called after the histogram is 
output — when the "CMD" prompt is displayed. If the user 
presses run/stop after the "CMD" prompt, the order 
statistics will automatically be displayed. 

The design of the program, especially the data entry 
loop, reflects the need for calculation speed. Code 
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redundancy exists at several points in order to reduce the 
need for extra flags, labels and goto statements which would 
slow execution during data entry. In spite of this need for 
speed, the summary totals needed for calculation of mean, 
standard deviation, skewness and kurtosis are accumulated 
during data entry. This is done so that these summary 
statistics are available with little or no wait following 
data entry. 

* 

A complete listing of the program registers and flags 
used by this program is shown at the end of the program 
listing. 

SAMPLE PROBLEM: 

In order to establish a training standard for an obsta- 
cle course, a division assistant 33 randomly selects 10 sol- 
diers and records the time it takes each to complete the 
course. The following times in minutes were recorded: 

2.1 2.4 2.2 2.7 2.5 

2.4 2.6 2.6 2.3 2.9 

Determine the summary statistics and cell frequencies 
necessary to plot a histogram of this data. 
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SOLUTION: 



1. First, set the size of the calculator's data memory 
large enough to retain the data values. This requires at 
least 16 registers plus 1 for each data point, or a total of 
26 in this example. Alternatively, the size of data memory 
may be set arbitrarily large, up to a maximum of 235 pro- 
vided the user has no other programs in the calculator he 
wishes to retain. For this example press: 

XEQ ALPHA SIZE ALPHA 26 

2. To call the program, determine the appropriate method of 
data entry and select the corresponding entry point. 

Press: 

XEQ ALPHA STAT1 ALPHA (quick entry) 

-or- 

XEQ ALPHA SI ALPHA (classroom use) 

3. The calculator will respond with the prompt "N?" asking 
for the number of data points. Press: 

10 R/S 

4. The calculator will respond with the prompt "XI?" asking 
for the first data point. Press: 

2.1 R/S 



49 



5. If you called the program via "si" the calculator will 
respond with "X1=2. 100" asking for verification that the 
first point is correct. If not correct enter the correct 
value, else press run/stop. 

6. The calculator will continue in the same way as steps 4 
and 5 for the remaining data points until all the data has 
been entered. If at any time you discover that you have 
made an error in data entry for any point, press: 

XEQ ALPHA SC ALPHA 

The calculator will respond with the prompt "POINT?" asking 
for the number of the point in error. For example, if point 
number 5 were in error, you would then press: 

5 R/S 

Assuming you had just input a 5 as the point in error, the 
calculator would then respond with the prompt "X5?" asking 
for the correct value of point 5. Respond with the correct 
value and press run/stop. The calculator will then go back 
to the place in the data entry sequence where it left off or 
it will go to the calculation of the summary statistics if 
data entry was previously completed. 
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7. 



When data entry has been completed, the calculator will 



respond with the mean of the data sample labeled as follows: 
XBAR=2. 470 



At this point, you have the option of reviewing the entire 
data set or continuing to calculate the remainder of the 
statistics. To review the entire data set, press: 



XEQ ALPHA SR ALPHA 



Note that if flag 21 is set on (press SP 21) , the calcula- 
tor will stop after each data point is displayed, permitting 
you to change any value simply by entering the new value and 
pressing run/stop. 



8. After the mean is displayed with the "XEAR" label, if 
you simply press the run/stop key, the calculator will cal- 
culate the following statistics with the label shown: After 

each press R/S. 



Di spl ay 



Meaning 



SSQD=0. 521 

SX=0. 24 1 
SKEW=0. 170 
KURT0= 2. 302 



Sum of Squared Deviations 
About the Mean 
Sample Standard Deviation 
Skewness 
Kurtosis 



9. At this point the calculator will automatically sort the 
data. This may take from several seconds to several minutes 
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depending on the number of points in the data set. After 
the data set has been sorted, the calculator will display 
the median as follows: 

MED=2 . 400 TO (Press H/S) 

. . 2.500 

Two data values are displayed because when the number of 
data points is even, the median is not unique, but rather 
spans an interval from the one point listed above to the 
other. Many users may wish to simply take the middle of 
this interval as the median, but any point is technically 
correct in the interval. When the number of data points is 
odd, the median is unique and only one value will be 
displayed by the calculator. 

10. After the median is displayed as described in step 9, 
the calculator will display the following statistics labeled 
as shown: 



Display 



Meaning 



HIN-2. 100 
MAX=2. 900 
RNG=0. 800 



Minimum Value 
Maximum value 
Range 



11. At this point the calculator will respond with ’’CELL?'* 
asking for the number of cells the user desires in the 
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histogram. If the number of cells is not significant at 
this point, the calculator will pick an appropriate number 
if the user simply presses run/stop. For this example, 
press: 

R/S 

12. Next the calculator responds with "WIDTH" asking for 
the width of the cells. Simply press run/stop if you do not 
wish to establish the width manually. Again, you may see 
the width the calculator will use by pressing the clear 
arrow key (Unless the width is an integer, you will also 
need to press FIX 3 to display the decimal properly if you 
wish to examine the width.) For this example, press: 

R/S 

13. The calculator will now display the cell frequency 
counts as an integer count followed by the next cell bound- 
ary. The leftmost cell boundary is set equal to the minimum 
value and is not explicitly output. If a data point falls 
exactly on a cell boundary, it is counted in the left cell. 
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For this example, the display will show: 



Dis£lay 

CNT=2 



Meaning 



xx=2. 260 
CNT=3 
XX=2 . 4 20 
CNT= 1 
xx=2 . 5 80 
CNT=3 
xx=2. 7 40 
CNT= 1 
xx=2. 900 



Two observations 

between 2. 1 (the minimum) 

and 2.26 (the cell boundary) 



Three observations 
between 2.26 (see above) 
and 2.42 (the next boundary) 



One observation 

between 2.42 (see above) 

and 2.58 (the next boundary) 



Three observations 
between 2.58 (see above) 
and 2.74 (the next boundary) 



One observation 
between 2.74 (see above) 
and 2.90 (the maximum) 



14. After the last cell boundary is displayed, the calcuia 
tor will display "CMD" asking the user for the next command 
Frequently, the user will wish to modify the histogram by 
changing the number of cells or the cell width. To recalcu 
late the histogram cell frequencies without re-entering the 
data press: 



If no further work with the histogram is desired, the user 
may view the order statistics simply by pressing run/stop. 



XEQ ALPHA AGAIN ALPHA 
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USER INSRUCTIONS : 



SINGLE VARIABLE STATISTICS 



;step 


EXPLANATION 


SEE 


PRESS 


RESULT 


: i 


SET SIZE 

(nnn= 16 + NUMBER 
OF DATA POINTS) 




XEQ 

"SIZE 

NNN 


UP TO 
nnr. 

= 235 


' 2 


CALL THE PROGRAM 
("STAT1 IS FOR 
REGULAR USE) 
("SI IS FOR 
CLASSROOM USE) 




XEQ 

"ST ATI 
-or- 
"S1 




; 3 


ENTER THE NUMBER 
DATA POINTS. 


N? 


input 

R/S 




' 4 


ENTER THE DATA 


XI?, X2? 
ETC. 


input 

R/S 






For mistakes or 
to review the data 
see last two steps 
below. 










WHEN VERIFY MODE 
IS SET ON (SET BY 
FLAG 05 ON) , AFTER 
EACH DATA POINT IS 
ENTERED, THE VALUE 
HILL 3E ECHOED BACK 
BY THE CALCULATOR. 


x1=xx 

etc. 


R/S 

-or- 

corr ect 
value 




; 5 


SUMMARY STATISTICS 
ARE CALCULATED WHEN 
ALL DATA HAS 3EEN 
ENTERED. 

STANDARD DEVIATION 

SKEWNESS 

KURTOSIS 


XBAR=xx 

SSQD=XX 

SX = XX 
SKEW =XX 
KURT=XX 


R/S- 

R/S 

R/S 

R/S 

R/S 


mean 

sum of 
sq dev 
from 
mean 


6 


CALCULATOR WILL 
AUTOMATICALLY SORT 
DATA POINTS. 


PRGM 




STANDBY 




AND THEN DISPLAY: 










MEDIAN 

(note rf N is even 
the median is not 
unique and an int- 
erval is displayed) 
MINIMUM 
MAXIMUM 
RANGE 


MED= XX 

MIN=XX 

MAX=XX 

RNG=XX 


R/S 

R/S 

R/S 

R/S 
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USSR IN 5RUCTI0NS 



SINGLE VARIABLE STATISTICS 



[STEP 


EXPLANATION 


SEE 


PRESS 


RESULT 


! 7 


USER OPTION TO ENTER 


CELL? 


R/S 






NUMBER OF HISTOGRAM 




-or- 






CELLS. NO INPUT IS 




INPUT N 






REQUIRED. 




R/S 




; 8 


USER OPTION TO ENTER 
WIDTH OF HISTOGRAM 
CELLS. (HAS PRECE- 
DENCE OVER NUMBER OF 
CELLS IF A WIDTH IS 
ENTERED. ) 


WIDTH? 


R/S 

-OR- 

INPUT 

R/S 






9 


CALCULATE HISTOGRAM 


CUT* II 


R/S 


CELL 




(OUTPUT DATA ABOUT 
EACH CELL FROM 
LEFT TO RIGHT.) 




FREQ 

COUNT 






XX = XX 


R/S 


UPPER 

X-VALUE 

LIMIT 






10 


ACCEPT NEXT COMMAND 


CMD 


ENTER 










NEXT 










CMD 




11 


RECALCULATE 




XEQ 






HISTOGRAM 




"AGAIN 




: 12 


EDIT AN INPUT VALUE 




XEQ 






AT ANY TIME PRIOR TO 
DATA SORT. 




"SC 








POINT? 


INPUT 


WILL 








POINT 


REMOVE 








NUMBR 


POINT 






X? 


INPUT 

CORRECT 










VALUE 






AFTER INPUT OF NEW 
VALUE CALCULATOR 
WILL RETURN TO DATA 
INPUT OR CALCULATION 
OF SUMMARY STATS AS 
APPROPRIATE. 
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USER INSBUCTIONS 



SINGLE VARI A3LE STATISTICS 



(STEP 


EXPLANATION 


SEE 


PRESS 


RESULT 


13 


REVIEW DATA POINTS 




XEQ 






(OR REVIEW ORDER 




"SR 






STATS AFTER SORT. ) 
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HP41C SOURCE CODE: 



SINGLE VARIABLE STATISTICS 



STAT1 



1 L3L "STAT1 

2 CF 05 

3 SF 26 

4 GTO "SS 



S 1 



RECOMMENDED ENTRY POINT 
SET VERIFY MODE OFF 
ENABLE AUDIO 



5 

6 
7 



LSL "SI 
SF 05 
CF 26 
8 SF 21 
( 



9 LSL "SS 

10 CF 29 

11 SREG 10 

12 CF 06 

13 CF 08 



14 L3L 00 

15 15 

16 STO 04 

17 STO 00 

18 RCL 15 

19 CIS 

20 "N? 

21 FC? 06 

22 PROMPT 

23 1E3 

24 / 



25 

26 



"SS 



ENTRY POINT FOR CLASSROOM USE 

SET VERIFY MODE ON 

DISABLE AUDIO TONES 

SET TO STOP DURING VERIFICATION 




27 STO 01 



ENTRY POINT FOR USSR SET OPTIONS 
NO DIGIT GROUPING 
ESTABLISH STATISTICAL REGISTERS 
USED BY DATA REVIEW FUNCTION 
‘USED BY DATA EDITING FUNCTION 



DATA ENTRY (F06-CLEAR) 
AND 

DATA REVIEW (F06-SET) 



(ESTABLISH INDIRECT ADDRESS BASE REG. 
(INITIALIZE DATA ENTRY POINTER 
^NUMBER DATA POINTS (LAST PROBLEM) 

(CLEAR MEANS DATA ENTRY, NOT REVIEW 



(SET UP LOOP COUNTER FOR DATA POINTS 
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HP41C SOURCE CODE 



SINGLE VARIABLE STATISTICS 



( 



28 


LBL 


01 




29 


ISG 


00 




30 


LBL 


02 




31 


RCL 


IND 


00 


32 


"X 






33 


FIX 


0 




34 


ARCL 


, 01 




35 


FIX 


3 




36 


AS TO 


i 03 




37 


FS? 


06 




38 


GTO 


03 




39 


»«-. = ? 






40 


TONE 


; 9 




41 


PROMPT 




42 


STO 


IND 


00 


43 


FC? 


05 




44 


GTO 


04 




45 


LBL 


03 




46 


CL A 






47 


ARCL 


, 03 




48 








49 


ARCL 


, IND 


00 


50 


CF 22 




51 


AVIEW 




52 


FC? 


22 




53 


GTO 


04 




54 


STO 


IND 


00 


55 


GTO 


03 




56 


LBL 


04 




57 


ST + 


10 




58 


XI 2 






59 


ST + 


1 1 




60 


LASTX 




61 


* 






62 


ST + 


12 




63 


LASTX 




64 


* 






65 


ST + 


13 




66 


FS? 


08 




67 


RTN 






68 


ISG 


01 




69 


GTO 


01 




70 


RCL 


01 





! 



DATA ENTRY LOOP 



(INCREMENT DATA STORAGE POINTER 
(RECALL DATA VALUE 



TEMP STORAGE FOR LABEL 

IS THIS REVIEW OF DATA PREV. ENTERED? 



(PROMPT USER FOR NEXT DATA VALUE 

(STORE THE DATA VALUE 

(NO VERIFICATION OF DATA DESIRED? 

(FOLLOWING IS THE VERIFICATION ROUTINE 

(RECALL THE LABEL 

( RECALL THE STORED DATA 
CLEAR DATA ENTRY FLAG 
WILL STOP FOR DATA ENTRY IF F2 1 SET 
WAS THERE NO DATA CHANGE DURING VIEW? 

(IF THERE WAS A NEW VALUE, THEN RECORD 
'IT AND GOBACK AND RE-VERIFY THE DATA. 
(FOLLOWING IS THE STATISTICAL ACCUM. 
'STORES SIGMA X 

(STORES SIGMA X -SQUARED 



(STORES SIGMA X-CUBED 



(STORES SIGMA X - FOURTH-PO W ER 
(IS THIS A DATA REVIEW? 

(IF DATA ENTRY, INCREMENT INPUT CNTR. 

(AT END OF DATA ENTRY, RECALL INPUT 
(COUNTER, WHICH IS A NUMBER EQUAL ONE 
(MORE THAN NUMBR POINTS 
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SINGLE VARIABLE STATISTICS 



( 



" SM 



CALCULATION OF SUMMARY STATS 



71 


LI 


!l 


"S 


72 


INT 




73 


1 






74 


— 






75 


STO 


15 


76 


MEAN 


77 


"XEAR 


78 


XEQ 


97 


79 


STO 


03 


80 


RCL 


1 1 


81 


RCL 


03 


82 


X| 2 




83 


RCL 


15 


84 


* 






85 


— 






86 


STO 


09 


87 


"SS CD 


88 


XEQ 


97 


89 


RCL 


15 


90 


1 






91 


— 






92 


/ 






93 


SQRT 


94 


"SX 




95 


XEQ 


97 


96 


RCL 


09 


97 


RCL 


15 


98 


/ 






99 


STO 


05 


100 


RCL 


12 


101 


RCL 


11 


102 


RCL 


03 


103 


* 






1 04 


3 






105 


* 






106 


- 






107 


RCL 


15 


108 








109 


RCL 


03 


1 10 


3 






111 


Y 1 X 




1 12 


2 






113 


* 






1 14 


+ 






115 


STO 


06 


116 


RCL 


05 


1 17 


1 . 


5 




1 18 


Y| X 




119 


/ 






120 


"SKEW 


121 


XEQ 


97 


122 


RCL 


13 


123 


RCL 


12 


124 


RCL 


03 



ENTRY ASSUMES X-REGISTER HAS A NUMBR 
1 MORE THAN NUMBER OF DATA POINTS. 



(STORES THE NUMBER OF DATA POINTS 



(CALL AN OUTPUT LA3ELING ROUTINE 
(TEMP STORE FOR XBAR 
’RECALL SIGMA X-SQUARED 
(RECALL XBAR 

(RECALL NUMBR POINTS 



(TEMP STORE FOR SUM OF SQUARED 
(DEVIATIONS ABOUT THE MEAN 

(NUMBER POINTS 

(CAN NOT USE SDEV FUNCTION BECAUSE OF 
(NON-STANDARD USE OF REGISTERS 12-14 

(STANDARD DEVIATION 

(SUM OF SQ DEVIATION ABOUT MEAN 
(NUMBER POINTS 

(SECOND MOMENT 
(SIGMA X-CUBED 
(SIGMA X-SQUARSD 
(XBAR 



(NUMBER POINTS 
(XBAR 



THIRD MOMENT 
SECOND MOMENT 



(OUTPUT THE SKEWNESS OF THE DATA 
(SIGMA X-FOURTH-POWER 
(SIGMA X-CUBED 
rx BAR 
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HP41C SOURCE CODE 
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125 


* 




126 


4 




127 


* 




128 


- 




129 


RCL 03 




130 


Xi 2 




131 


RCL 11 




132 


* 




133 


6 




134 


* 




135 


+ 




136 


RCL 15 




137 


/ 




138 


RCL 03 




139 


4 




140 


Y | X 




141 


3 




142 


* 




143 


— 




144 


STO 07 




145 


RCL 05 




146 


X| 2 




147 


/ 




148 


"K URT= 




149 


XEQ 97 
(FS? 09 
l RTN 




1 50 


XEQ 98 




151 


CF 00 




152 


RCL 15 




153 


2 




154 


/ 




155 


FRC 




156 


X=0? 




157 


SF 00 




158 


LASIX 




159 


.5 




160 


+ 




161 


RCL 04 




162 


+ 




163 


"MEE = 




164 


ARCL IND 


X 


165 


FC? 00 




166 


GTO 05 




167 


TO 




168 


PROMPT 




169 


1 




170 


+ 




171 


it it 

« « 




172 


ARCL IND 


X 


173 


L3L 05 




174 


FROMPT 





(XBAR 

(SIGMA X-SQUARED 



(NUMBER POINTS 
(XBAR 



(FOURTH MOMENT 
(SECOND MOMENT 



(OUTPUT THE KURTOSIS 
(SHORT FORM WOULD NOT COMPUTE STATS 
(WHICH REQUIRE SORTED DATA 
(CALL A DATA SORTING ROUTINE 
(INITIALIZE TEMP FLAG USED TO CHECK 
(EVEN OR ODD NUMBER OF DATA POINTS 



(WAS IT AN EVEN NUMBER OF POINTS? 
(IF WAS EVEN NUMBER, SET FLAG. 

(COMPUTING ADDRESS OF MEDIAN 

(ADDRESS BASE REGISTER 

(X-REG NOW HAS ADDRESS OF MEDIAN 



(EVEN NUMBER POINTS IMPLIES THE MEDIAN 
(NOT UNIQUE, BUT SPANS AN INTERVAL 

(DISPLAY THE LEFT BOUNDARY OF MEDIAN 

(X-REG POINTS TO RIGHT BOUND OF MEDIAN 
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HP41C SOURCE CODE 



SINGLE VARIABLE STATISTICS 



( 



"AGAIN 



DISPLAY HISTOGRAM 



175 


LBL 


"AGAIN 


176 


RCL 


04 




177 


1 






178 


+ 






179 


RCL 


IND 


X 


180 


"MIN 




181 


XEQ 


97 




182 


STO 


09 




183 


RCL 


04 




184 


RCL 


15 




185 


+ 






186 


RCL 


IND 


X 


187 


"MAX 




188 


XEQ 


97 




189 


RCL 


Z 




190 


— 






191 


STO 


03 




192 


"RNG 




193 


XEQ 


97 




1 94 


CF 00 




195 


RCL 


15 




196 


RCL 


04 




197 


+ 






198 


1E3 






199 


/ 






200 


RCL 


04 




201 


+ 






202 


1 






203 


+ 






204 


STO 


01 




205 


RCL 


15 




206 


LN 






207 


2 






208 


* 






209 


FIX 


0 




210 


RND 






211 


"CELL? 




212 


PROMPT 




213 


RCL 


08 




214 


XOY 




215 


/ 






216 


"WIDTH 




217 


PROMPT 




218 


STO 


08 




219 


LBL 


06 




220 


RCL 


08 




221 


ST + 


09 




222 


0 






223 


STO 


02 




224 


LBL 


07 




225 


RCL 


IND 


01 


226 


RCL 


09 




227 


FIX 


3 




228 


RND 







(ADDRESS BASE REGISTER 



(RECALL THE FIRST ORDER STAT 

(CALL AN OUTPUT LABELING ROUTINE 
(HOLDS STARTING (LEFTMOST) X BOUNDARY 
(ADDRESS BASE REGISTER 
(NUMBER OF DATA POINTS 

(RECALL THE N-TH ORDER STATISTIC 

(DISPLAY THE MAX VALUE OBSERVED 
(MIN 

(TEMP STORE FOR THE RANGE 
(DISPLAY THE RANGE 

(INITIALIZE TEMP FLAG TO MARK LAST BAR 
(NUMBER DATA POINTS 
(ADDRESS BASE REGISTER 



(COMPUTING INDEX LOOP COUNTER 
(ADDRESS BASE REGISTER 



R 01 SET TO ADDRESS AND LOOP THRU DATA 
NUMBER POINTS 



(DEFAULT NUMBER OF BARS IS 2*LN (N) 

(VALUE IS ROUNDED NOT TRUNCATED 

(USER HAS OPTION TO CHNGE NUMBR CELLS 
(RANGE 



(USER HAS OPTION TO CHANGE CELL WIDTH 
(NOW HOLDS CELL WIDTH NOT RANGE 

(CELL WIDTH 

(UPPER LIMIT OF CURRENT CELL COUNTED 

(INITIALIZE CELL COUNTER 

(NEXT DATA POINT 
(CELL UPPER LIMIT 
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HP 4 1 C SOURCE CODE 



SINGLE VARIABLE STATISTICS 



229 


X< Y? 




230 


GTO 


08 


231 


1 




232 


ST+ 


02 


233 


IS G 


01 


234 


GTO 


07 


235 


SF 00 


236 


LBL 


08 


237 


"CNT 




238 


RCL 


02 


239 


FIX 


0 


240 


XEQ 


97 


241 


FIX 


3 


242 


"XX 




243 


RCL 


09 


244 


XEQ 


97 


245 


FC?C 


00 


246 


GTO 


06 


247 


"CMD 


i 


248 


AVIEW 


249 


RTN 





(DATA POINT LESS THAN UPPER LIMIT 

(INCREMENT THE CELL COUNTER 
(PREPARE TO LOOK AT NEXT DATA POINT 

(SET FLAG FOR OUTPUT OF LAST BAR 



(OUTPUT THE CELL FREQUENCY COUNT 



OUTPUT CELL BOUNDARY — LOWER LIMIT 
IS THIS THE LAST BAR ? 
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HP41C SOURCE CODE 



SINGLE VARIABLE STATISTICS 



( ( 

( "SR 



250 


LBL 


"SR 


251 


SF i 


06 


252 


GTO 


00 








253 


LI 


JL 


"SC 


254 


SF 


08 


255 


CF 1 


06 


256 


RCL 


00 


257 


STO 


05 


258 


RCL 


01 


259 


STO 


06 


260 


"POINT? 


261 


FROMPT 


262 


STO 


01 


263 


RCL 


04 


264 


+ 






265 


STO 


00 


266 


RCL 


IND 00 


267 


S3 


1 _ 


10 


268 


XI 


'2 




269 


S3 


1 _ 


1 1 


270 


LASTX 


271 


* 






272 


ST- 


12 


273 


LASTX 


274 








275 


ST- 


13 


276 


XEQ 


02 


277 


CF ' 


08 


278 


RCL 


05 


279 


STO 


00 


280 


RCL 


06 


281 


STO 


01 


282 


1 






283 


— 






284 


IS G 


X 


285 


GTO 


02 


286 


GTO 


"SM 



REVIEW THE DATA 



(SETS MODE FOR REVIEW NOT QUERY 



EDIT THE DATA 



SET TO EDIT MODE 

SET TO QUERY FOR CORRECT VALUE 

CURRENT INPUT ADDRESS POINTER 

SAVE TO ENABLE RETURN TO DATA ENTRY 

CURRENT INPUT INDEX LOOP COUNTER 

SAVE TO ENABLE RETURN TO DATA ENTRY 



(ESTABLISH PSEUDO-INDEX COUNTER 
(ADDRESS BASE REGISTER 

(COMPUTED ADDRESS OF DATA TO 3E EDITED 
(RECALL THE OLD VALUE 
(CORRECT SIGMA X 

(CORRECT SIGMA X-SQUARED 



(CORRECT SIGMA X -CUBED 



(CORRECT SIGMA X-FOURTH-PO WER 
(CALL DATA ENTRY AS A SUBROUTINE 
(FOLLOWING RESTORES DATA ENTRY 
(INPUT ADDRESS REGISTER VALUE 

(INPUT LOOP COUNTER 



(TEST TO SEE IF ALL DATA ALREADY INPUT 
'IF NOT, BRANCH TO THE INPUT LOOP 
(IF YES, RECOMPUTE THE SUMMARY STATS 
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HP41C SOUBC E CODE 



SINGLE 7ARIA3LE STATISTICS 



(( 



97 OUTPUT LABELING ROUTINE 



287 L3L 97 

288 

289 ARCL X 

290 PROMPT 

291 RTN 



6.5 



HP41C SOURCE CODE: 



SINGLE VARIABLE STATISTICS 



( 



98 



SHELL SORT 



292 


LE 


II 


98 




293 


RCL 


15 




2 94 


STO 


01 




295 


L3L 


09 




296 


RCL 


01 




297 


2 








298 


/ 








299 


INT 






300 


STO 


01 




301 


X= 


= 0? 






302 


RTN 






303 


1 








304 


STO 


02 




305 


L3L 


10 




306 


STO 


03 




307 


LBL 


1 1 




308 


RCL 


01 




309 


♦ 








310 


RCL 


04 




311 


+ 








312 


RCL 


IND 


X 


313 


RCL 


03 




3 14 


RCL 


04 




315 


+ 








316 


XOY 




317 


RCL 


IND 


Y 


318 


X< = Y? 




319 


GTO 


12 




320 


STO 


IND 


T 


321 


XOY 




322 


STO 


IND 


Z 


323 


RCL 


03 




324 


RCL 


01 




325 


- 








326 


STO 


03 




327 


X>0? 




328 


GTO 


11 




329 


LBL 


12 




330 


RCL 


15 




331 


RCL 


01 




332 


- 








333 


RCL 


02 




334 


1 








335 


+ 








336 


STO 


02 




337 


X< = ’i 






338 


GTO 


io 




339 


GTO 


09 





(RECALL NUMBER OF DATA POINTS 
(DEFINE A = "MIDPOINT" 

(RECALL MIDPOINT 



(A = I NT (A/2 ) 

(TEST TO SEE IF LIST SORTED 



B = 



1 -- RESET LEFT SHELL BOUNDARY 

STACK TABLE FOLLOWS: 



C = B 
C 
A 

D=C+ A 
BASE 
ADDR I 

BASE 
ADDR ( 

x Sci 



B 
B 
C 
B 
D 
B 

ADDR 

X< D > 

X (D) 
ADDR 
X (D) 



B 

ADDR D 

Hffi » 

ADDR D 
ADDR C 



,B 

, ADDR D 



, A DDR D 



[FOLLOWING INTERCHANGES X (C) AND X(D) 

ADDR C , ADDR 




ADDR 

ADDR 



C 

C 

C 



, ADDR 
, ADDR 



D 

D 

D 



N 

A N 

E=N- A 

B E 

1 3 

B + 1 E 

8 = 8+1 E 



E 
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HP41C SOURCE CODE 



SINGLE VARIABLE STATISTICS 



340 END 



THIS PROGRAM USES THE FOLLOWING REGISTERS: 

R00 — INPUT DATA ADDRESS POINTER 
RO 1 — LOOP INDEX COUNTER 

(USED BY DATA ENTRY AND SORT ROUTINES) 
R02 — TEMP REGISTER 

(CELL FREQUENCY COUNT IN HISTO RTN) 
j AND SHELL BOUNDARY IN SORT ROUTINE) 

RO 3 — TEMP REGISTER 

i INPUT LABEL IN DATA INPUT ROUTINE) 

XBAR IN SUMMARY STAT ROUTINE) 

AND POINTER IN SORT ROUTINE) 

NDIRECT ADDRESS BASE 

RO 5 — SECOND MOMENT (POPULATION VARIANCE) 

R06 — THIRD MOMENT 
R07 — FOURTH MOMENT 
RO 8 — HISTOGRAM CELL WIDTH 
R09 — TEMP REGISTER 

(HOLDS SUM OF SQ DEVIATION ABOUT MEAN) 
(AND HISTOGRAM CELL UPPER LIMIT) 

RIO — SUM OF X VALUES 
R 1 1 — SUM OF X-SQUARED VALUES 
R 1 2 — SUM OF X-CUBED VALUES 
R 1 3 — SUM OF X RAISED TO THE FOURTH POWER 
R 1 4 — NOT USED BUT SET TO ZERO BY CLRS 
R 1 5 — NUMBER DATA POINTS (SET BY &+) 
R16....R228 RAW DATA POINTS 

— IN NATURAL SEQUENCE BEFORE SORT 
— AS ORDER STATISTICS AFTER SORT 

THIS PROGRAM USES THE FOLLOWING FLAGS: 

FOO — TEMP FLAG (USED IN EDIT AND HISTO RTNS) 
FO 5 — VERIFY MODE (EVERY DATA POINT ECHOED) 
F06 — INDICATES REVIEW OF DATA NOT QUERY MODE 
F08 — INDICATES EDITING A DATA POINT 
?2 1 — PRINTER ENABLE (STOPS CALCULATOR 
DURING A VIEW INSTRUCTION) 

F26 — AUDIO ENABLE 
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SINGLE VARIABLE STATISTICS 



13 



14 

■ 
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in 

iiiiniiiiiiiiiiiiiiiiDiiiiiiiiniiiii 




llillJillliliailllltllillllffl 




iffiiiiHiniiiiiiifiiniiiniim 






jllll llllSMailll ij HI U li l iBII 
C lII I I I IBIlBBill l lBIMIIIillllilBill 



aaijiiMiiiiMJimu 




mam i 





53 



SINGLE V ft R I ft 8 l E STftTISTICS 




llllllllllllllllilllllllMlllllllimffliyillllMtllBIBIIlM 



27 








20 












— — 





70 




SINGLE VARIABLE STATISTICS 

37 



38 




in 



39 



40 



41 



4? 



43 



44 






























■ 




















i 
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APPENDIX 3 

LINEAR PROGRAMMING EXAMPLE 

INTRODUCTION: 

Linear programming is a n operations research technique 
normally associated with computerized data bases and the 
largest computers. Because of the complexity of the com- 
puter programs for linear programming and the large amount 
of data associated with most real world problems, calcula- 
tors have not been widely used for this application. With 
the increased memory capacity of the HP41CV, however, it is 
now possible to offer a calculator program which can solve 
interesting small scale linear programs. Of value primarily 
as an educational aid, this program will also be able to 
solve many small scale problems found at Army division, bri- 
gade and battalion level. The text by Hillier and Lieberman 
[Ref. 8: pp. 16-66] is recommended as an introduction to 

the theory of linear programming as used by the program 
given in this appendix. Use of the program requires the 
user to formulate the linear programming problem; set up a 
Simplex tableau in standard form including adding slack, 
surplus and artificial variables as required; and interpret 
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rhe final tableau including the calculation of the values 
associated with the variables in the final basis. Using the 
tableau form of the Simplex algorithm, the calculator per- 
forms both phase I (to obtain a feasible solution) and phase 
II (to obtain an optimal solution) to solve the linear pro- 
gramming problem. The calculator automatically determines 
the pivot column and pivot row for each pivot step. Infeas- 
ible and unbounded problems are automatically identified for 
the user by the program. There is no explicit handling of 
variables with upper bounds. 

PROGRAM DESCRIPTION: 

The program is written as a series of subroutines, each 
of which performs a major step in the Simplex algorithm. To 
provide clarity to the user, alphabetic labels have been 
retained to identify the subroutines in lieu of faster and 
more memory efficient numeric labels. The alphabetic labels 
have not been retained for use as program entry points and 
may be changed to numeric labels at the option of the user. 
The program has two entry points, "LP" for running a new 
problems and "ALP” for reviewing data previously entered. 

Subroutine "FINDQ" determines the pivot column by 
selecting the variable in the objective function with the 
most negative "price." If "FIN DQ" discovers at least one 
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negative value in the objective function, then the tableau 
column number associated with the most negative value will 
be stored in register 05. Upon return from "FINDQ , " the 
main routine tests register 05 to see if it contains a non- 
zero entry. If the entry is zero, it means that no further 
pivots will improve the value of the objective function, and 
the Simplex algorithm halts. If the program was in phase I 
(flag 1 1 clear) and the value of the phase I objective func- 
tion is reduced to zero, then a feasible solution has been 
found and the program will automatically proceed to phase II 
to discover an optimal solution. 

Subroutine "FINDP" determines which variable will leave 
the current basis by performing a minimum positive ratio 
test along the pivot column. In this way, the pivot row is 
determined. The row number of the pivot row is stored in 
register 06. Upon return from subroutine "FIN DP , " the main 
routine tests register 06 to see if it contains a non-zero 
entry. If the entry is zero, it means that the problem is 
unbounded and the Simplex algorithm halts. Such an 
unbounded condition is most likely caused by an error in the 
problem formulation. 

Having determined the pivot column and the pivot row, 
subroutine "PIVOT" performs the actual Simplex pivot 
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operation. To speed calculation register 00 is used as a 
temporary register to hold the reciprocal of the pivot 
element. Note that the pivot row is handled separately from 
the other rows in the tableau. Flag 0U is used to provide 
the option of stopping calculation after every pivot. When 
this flag is set, the program will halt to allow the user to 
review the status of the tableau with the "ALP" function. 

Subroutine "CHECK" has two primary functions. First, it 
is used to verify than the designated basic variables are in 
row elimination form prior to the start of the Simplex 
algorithm. This means that the basic variable must have a 
coefficient of 1 in the row in which it is basic and zero's 
in all other rows. The second function of check is to pre- 
pare the objective function for phase I, if the initial 
basis contains artificial variables as indicated by one or 
more minus signs in the "JB" vector. 

Three subroutines are used to query the user for input 
data. Subroutine "READMN" queries the user for the number 
of constraints and decision variables in the problem and 
verifies the calculator memory is set to contain all the 
data necessary to solve the problem. Subroutine "READJ3" 
queries the user for a column vector of pointers which indi- 
cate which variable is currently basic in each row. When 
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entering this vector of pointers, the user indicates 
artificial variables with a minus sign. Subroutine "READA" 
gueries the user for the values in the initial Simplex 
tableau including the slack and surplus variables and the 
right hand side and objective function. 

Several other service routines also are provided in this 
program. Memory size verification is done by subroutine 
"SIZE, 1 ’ which is called from within "READMN." Subroutine 
"IN” is used to query the user for data entry and is called 
by all of the data input routines. Subroutine ”NXT n ini- 
tializes registers which contain frequently used quantities 
such as the the total size of the tableau for phase I atid 
phase II. Subroutine "INIT" clears the calculator memory 
and sets flags and program constants appropriately for input 
of a new problem. Subroutine "SETL'’ establishes the loop 
counters used repeatedly within almost every other subrou- 
tine. Subroutine " ERR1 " displays an appropriate error 
message if a data entry error is detected. 

SAMPLE PROBLEM: 

A division assistant G4 is planning an ammunition upload 
plan. There are four types of tank munitions to consider. 
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incl tiding: 



A = Discarding Sabot Hounds 
3 = High Explosive Anti-Tank Hounds 
C = Phosphorous Munitions 
D = Machine Gun Ammunition 



Based on the Commander's guidance the assistant G4 is to 



consider the sabot rounds as twice as important as the HEAT 



rounds, which in turn are themselves twice as important as a 



unit amount of phosphorous munitions and machine gun 



ammunition. His mission then, is to maximize: 



Z=4A+2B+C+D 



He is, however, constrained by the following factors: 



1. There can be no more than 30 units of both sabot 
and HEAT munitions combined. 

2. There can be no more than 50 units of all types of 
ammunition combined. 

3. There must be at least 30 units of HEAT and 
phosphorous munitions combined. 

4. There must be at least 5 units of machine gun 
ammunition. 

These constraints may be expressed as: 

A + B <30 

A + B+C+D < 50 

B + C >30 

D > 5 



Based 01 * the Commander's guidance and the constraints listed 



above, formulate an optimum load plan. Fractional units are 



permitted. 
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SOLUTION: 



1. Before beginning with the calculator, the first step is 
to layout the tableau in standard form. This step and the 
last step of interpreting the final tableau require working 
knowledge of linear programming as explained in Hillier and 
Liberman [Ref. 8: pp. 16-66]. The standard form of the 

tableau is: 



1 2 3 

A 3 C 

1 1 

1 1 1 

1 1 

-ztj r-2 -T~ 



4 5 6 

D HI H2 

1 

-1 

J ’I 



7 8 9 

SI S2 A 1 

1 

1 

1 



10 11 

A2 RHS 

30 

50 

30 

1 5 

u — 



In this tableau, HI and H2 are surplus variables; SI and S2 
are slack variables; and A1 and A2 are artificial variables. 



2. The first step with the calculator is to set the size of 
the calculator's data memory. This program requires 20 reg- 
isters for temporary storage, 1 register for each tableau 
element, and 1 register for each row to hold the pointer to 
the basic variable for that row. Thus, if M is the number 
of constraints and N is the number of variables including 
slack, surplus and artificial variables, then the total data 
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storage requirement is: 



storage required = 21 ♦ M ♦ ( (N ♦ 1) x (H + 2) ) 

As mentioned in the program description, the "SIZE" subrou- 
tine will automatically verify that the user has allocated 
enough data storage to solve the problem. The length of the 
program is such that 177 data storage registers is the maxi- 
mum number of data storage registers that can be allocated. 
Thus, linear programs with 7 constraints and 15 decision 
variables can be solved with this program. For this 
example, press: 

■* XEQ ALPHA SIZE ALPHA 175 

3. Call for execution of the program with a new data set. 
Press: 



XEQ ALPHA LP ALPHA 

4. The calculator will respond with the prompt "MOM ROWS?” 
asking for the number of constraints in the linear program 
formulation. In this example, press: 

4 R/S 

5. The calculator will respond with the prompt "NUM COLS?” 
asking for the number of variables in the problem. The user 
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must count the number of slack, surplus and artificial 
variables in this total. In this example, press: 

10 R/S 

6. The calculator will respond with the prompt "BASIC 1 ?" 
asking for the variable number of the variable which is 
basic in the first row. One of the major features of this 
program is that the basic variables need not be in the 
rightmost positions in the tableau. Thus, if a tableau were 
given in which some pivots had already been performed, the 
program could resume operation immediately. In this 
example, press: 

7 R/S 

In a similar fashion, the calculator will then query the 

user for the variable number of the variables which are 

basic in the remaining rows. 

For this example: 

See Respond With 

Basic 2 ? 8 R/S 

Basic 3 ? 9 CHS R/S 

Basic 4 ? 10 CHS R/S 

Notice that because the basic variables in rows three and 
four are artificial variables, the variable number is 
entered as a negative number. This signals the calculator 
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that these variables must be driven from the basis in order 



to reach an initial feasible solution. 

7. Next, the calculator will respond with "T1,1?" asking 
for the first element in the tableau. The user should enter 
the numbers in the tableau using the digit entry keys and 
pressing run/stop after every entry. Notice that the right 
hand side and the objective function will be entered with 
the appropriate index in the tableau as shown in step 1 
above. The user must insure that the objective function is 
in standard form with the appropriate sign for each 
coefficient — in this example each coefficient is negative. 

8. After the last element in the tableau has been entered, 
the calculator will begin to automatically perform the Sim- 
plex algorithm. If the user wishes to stop the calculator 
after every pivot, he may at any time press: 

R/S SF 04 R/S 

If this flag is set, the calculator will stop and display 
the pivot number after every pivot is completed. 

9. When the Simplex algorithm can no longer improve the 
objective function, the calculator will stop and display the 
value of the objective function. In this example, the 
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calculator will stop after approximately three minutes and 



display : 



VALtJE=1 10.000 



10. At this point, the user must use entry point "ALP" to 
determine the status of the final tableau. For this 
example, press: 

XEQ ALPHA ALP ALPHA 



Then by sequentially pressing the run/stop and clear arrow 
keys, the basic variables and tableau entries will be 
displayed. For example, in this problem; 

See P r es s See Meaning 

BASIC 1? C LX 2 Variable 2 is basic 

in the first row. 

BASIC 2? CLX 1 Variable 1 is basic 

in the second row. 



etc. 



Then for the elements of the tableau: 

See Press See Meaning 

T 1 , 1 ? CLX 0.000 Tableau entry 

T 1 , 2? CLX 1 .000 Tableau entry 

11. After the calculator is finished, it remains for the 

user to interpret the final tableau. Again, the reference 

by Hillier and Lieberman [Ref. 8: pp. 16-66] is of primary 

value. In particular, the user must be able to determine 

the value of the final decision variables based upon what 
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variables are in the basis, and vnat the final "right hand 
side" values are for each row. For this example, the final 
tableau is: 



3 4 5 

C D HI 

-1 

1 

1 

1 

— 0 TJ U 3 IT 

Thus, the solution may be 
munitions A,B and C and 5 



1 2 

A B 

1 

1 



6 


7 


8 


9 


10 


1 1 


H2 


SI 


S2 


A 1 


A 2 


RHS 


- 1 


1 


-1 


1 


1 


15 


1 




1 


>1 


-1 


15 


1 




1 




-1 


15 


-1 








1 


5 


7 


~T 


3 


-7 


~~2 


TTU 


.nter pr 


s ts d 


as 15 


uni 


ts each 


for 



units for munition D. 
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ISER INDUCTIONS: 



LINEAR PROGRAMMING 



( 

(STEP 

( 


EX PLA N AT ICN 


SEE 


PRESS 


RESULT 


( 

( 1 
( 

( 

( 

( 


SET SIZE 
( NNN= 

2 1 + M+ ( N + 1 ) (M + 2 ) 
WHERE M=N JM ROWS 
AND N-NJM COLS 




XEQ 

■'SIZE 

NNN 


UP TO 
N IN 

= 177 


- 

( 

( 2 
( 

( 


CALL THE PROGRAM 




XEO 

"LP 




( 

( 2 
( 

( 


ENTER THE NUMBER 
CF CONSTRAINTS 


NUM ROW 
? 


INPUT M 
R/ S 




( 

( 4 
( 

( 

( 


ENTER THE NUMBER 
CF VARIABLES 
( INCLUDE SLACKS , 
SURPLUS C ARTIF. 


NUM COL 
? 


INPUT N 
R/S 




( 

( 5 
( 

( 

( 


ENTER CURRENT 
BASIC VARIABLE 
NUMBERS BY ROW 


BASIC I 
? 

ETC . 


INPUT 
VAR V 
R/S 




( 

( 6 
( 

( 

( 

( 

( 

( 

( 


ENTER TABLEAU 
VALLES. 

FOR MISTAKES OR 
TO REVlE'w THE DATA 
SEE LAST STEP 
BELCN. 


TI, 1? 

ETC . 


INPUT 

R/S 

R/S 
SF OA 
R/S 


— 


( 

( 7 
( 

( 

( 

( 


TO FCRCE THE CAL- 
CULATOR TO STOP 
AFTER EACH PIVOT. 




< 

( 6 
( 

( 

( 


SIMPLEX COMPLETED: 
OPTIMAL SOLUTION 
FOUND. 


VALUE= 
XX. XXX 




F INAL 
OBJ. 

FUNCTION 

VALUE 


( 

( S 
( 

( 

( 


SIMPLEX COMPLETED: 
PROBLEM IS 
INFEASIBLE. 


INFEAS 
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USER INSFUCT ICNS: 



LINEAR PROGRAMMING 



( 

(STEP 
{ 


EXPLANAT ICN 


SEE 


PRESS 


RESULT 


( 

( IG 


SIMPLEX COMPLETED: 


UNBOUND 






I 


PROBLEM IS 








( 

( 

{ 


UNBOUNDED. 








( 

( 11 


TO REV I Eh VALUES 




XEQ 




{ 


IN TABLEAU AT ANY 




"ALP 




( 


TIME, INCLUCING 








{ 

( 


FINAL TABLEAU. 








( 

( 

( 

( 

( 

{ 

( 


AS PRGMPTS APPEAR, 
DATA CAN BE CHANGED 
BY ENTERING NEW 
VALLE. 








WHAT IS CURRENTLY 


BASIC 1 


CLX 


PROMPT 


( 


BASIC? 


-> 




* ILL 


( 




ETC . 




VANISH 


{ 








LEAV ING 


( 

( 

{ 








CATA 


WHAT ARE VALUES 


T1 , 1 ? 


CLX 


PROMPT 


( 


IN TABLEAU? 


ETC . 




WILL 


I 








VANISH 


( 








LEAVING 


( 


OBTAIN VALUES 






DATA 


( 


CF FINAL SOLU- 








( 


TION FROM 








( 


KNOWING WHICH 








( 


VARS ARE BASIC 








( 


ANC VALUE GF 








{ 


RIGHT-HAND- 








( 


SICE FRCM THE 








( 

I 


TABLEAU. 
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HP41C SOURCE CODE 



LINEAR PROGRAMMING 





| (- 




1 


LBL 


"L? 


2 


XEQ 


"INIT 


3 


LBL 


"ALP 


4 


XEQ 


"READMN 


5 


XEC 


"READJB 


6 


XEQ 


"REACA 


7 


XEQ 


"CHECK 


8 


LBL 


15 


9 


XEQ 


"FINDQ 


10 


RCL 


05 


1 1 


X# 0? 




12 


GTO 


35 


13 


FS? 


1 1 


14 


GTO 


25 


15 


RCL 


04 


16 


RCL 


10 


17 


RCL 


12 


18 


* 




19 


+ 




20 


RCL 


07 


21 


RCL 


IND Y 


22 


A3S 




23 


X< Y? 


24 


GTO 


20 


25 


"INFSAS 


26 


AVIEW 


27 


STOP 


28 


LBL 


20 


29 


SF 1 1 


30 


RCL 


04 


31 


RCL 


08 


32 


RCL 


12 


33 


* 




34 


+ 




35 


STO 


14 


36 


GTO 


15 


37 


LBL 


25 


38 


RCL 


04 


39 


RCL 


09 


40 


RCL 


12 


41 


* 




42 


+ 




43 


RCL 


IND X 


44 


STO 


00 


45 


"VALUE= 


46 


ARC I 


, 00 


47 


AVIEW 


48 


STOP 


49 


LBL 


35 


50 


XEQ 


"FINDP 


51 


RCL 


06 


52 


X#0? 


53 


GTO 


40 


54 


"UNBOUND 
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HP41C SOURCI CODE 



LINEAR PROGRAMMING 



55 ARCL 05 

56 STOP 

57 L3L 40 

58 XEQ "PIVOT 

59 GTO 15 
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HP 4 1 C SOURCE CODE 



LINEAR PROGRAMMING 



( 



60 


LBL 


61 


SF 0 


62 


0 


63 


STO 


64 


STO 


65 


1 


66 


RCL 


67 


XEQ 


68 


LBL 


69 


RCL 


70 


RCL 


71 


+ 


72 


STO 


73 


RCL 


74 


RCL 


75 


— 


76 


RCL 


77 


CHS 


78 


X< Y? 


79 


GTO 


80 


FC? 


81 


GTO 


82 


RCL 


83 


RCL 


84 


+ 


85 


RCL 


86 


ABS 


87 


RCL 


88 


X< = Y 


89 


GTO 


90 


LBL 


91 


RCL 


92 


STO 


93 


RCL 


94 


INT 


95 


STO 


96 


LBL 


97 


ISG 


98 


GTO 


99 


CF 0 


1 00 


BIN 



FINDQ 



"FINDQ 

03 

05 

11 

"SETL 

31 

14 

01 

00 

IND X 

03 

07 

38 

11 

37 

15 

01 

IND X 

07 

38 

37 

IND 00 

03 

01 

05 

38 

01 

31 

1 



88 



HP 4 1 C SOURCE CODE: LINEAR PROGRAMMING 

; ( 

FINDP 



101 


LBL 


"FINDP 


102 


SF 


02 




103 


0 






104 


STO 


06 




105 


1E20 




106 


STO 


03 




107 


1 






108 


RCL 


08 




1 09 


XSQ 


"SETL 


1 10 


LBL 


41 




111 


RCL 


01 




112 


1 






113 


— 






114 


RCL 


12 




1 15 


* 






1 16 


RCL 


04 




1 17 


+ 






1 18 


STO 


00 




1 19 


RCL 


05 




120 


♦ 






121 


RCL 


IND 


X 


122 


STO 


02 




123 


RCL 


07 




124 


X>Y 


•> 

m 




125 


GTO 


48 




126 


RCL 


00 




127 


RCL 


12 




128 


+ 






129 


RCL 


IND 


X 


130 


RCL 


02 




131 


/ 






132 


STO 


00 




133 


RCL 


03 




134 


- 






135 


RCL 


07 




136 


CHS 






137 


X< Y 


■> 

• 




138 


GTO 


48 




139 


LBL 


47 




140 


RCL 


00 




141 


STO 


03 




142 


RCL 


01 




143 


INT 






144 


STO 


06 




145 


LBL 


48 




146 


IS G 


01 




147 


GTO 


41 




148 


CF 


02 




149 


RTN 
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HP 4 1 C SOURCE CODE 



LINEAR PROGRAMMING 



( 

PIVOT 



150 


L3L 


"PIVOT 


151 


SF 03 




152 


RCL 


06 




153 


1 






154 


— 






155 


RCL 


12 




156 


* 






157 


RCL 


04 




158 


+ 






159 


STO 


03 




160 


RCL 


05 




161 


+ 






162 


RCL 


IND 


X 


163 


1/X 






164 


STO 


00 




165 


1 






166 


RCL 


12 




167 


XEQ 


"SETL 


168 


L3L 


51 




169 


RCL 


03 




170 


RCL 


01 




171 


+ 






172 


RCL 


00 




173 


ST* 


IND 


Y 


174 


ISG 


01 




175 


GTO 


51 




176 


1 






177 


RCL 


09 




178 


FS? 


11 




179 


0 






180 


FC? 


1 1 




181 


1 






182 


+ 






183 


XEQ 


"SETL 


1 84 


L3L 


52 




185 


RCL 


06 




1 86 


RCL 


01 




187 


INT 






188 


X=Y? 




189 


GTO 


59 




190 


1 






191 


— 






192 


RCL 


12 




193 


* 






1 94 


RCL 


04 




195 


+ 






196 


STO 


16 




1 97 


RCL 


05 




198 


+ 






199 


RCL 


IND 


X 


200 


STO 


00 




201 


2 






202 


RCL 


12 




203 


XEQ 


"SETL 
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HP41C SOURCE CODE 



LINEAR PROGRAMMING 



204 


L3L 


53 


205 


CP 07 


206 


RCL 


05 


207 


RCL 


02 


208 


INT 




209 


X=Y? 




210 


SF 07 


211 


RCL 


03 


212 


♦ 




213 


RCL 


IND 


214 


RCL 


00 


215 


* 




216 


RCL 


16 


217 


RCL 


02 


218 


+ 




2 19 


XOY 


220 


ST- 


IND 


221 


FC? 


07 


222 


GTO 


54 


223 


0 




2 24 


STO 


IND 


225 


LBL 


54 


226 


ISG 


02 


227 


GTO 


53 


228 


LBL 


59 


229 


ISG 


01 


230 


GTO 


52 


231 


RCL 


13 


232 


RCL 


06 


233 


+ 




234 


RCL 


05 


235 


STO 


IND 


236 


1 




237 


ST + 


17 


238 


CF 03 


239 


TONE 


: 7 


240 


FC? 


04 


241 


RTN 




242 


"PI VOT 


243 


ARCI 


, 17 


244 


PROMPT 


245 


RTN 
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HP41C SOURCE CODE 



LINEAR PROGRAMMING 



( 

READMN 



246 


LBL 


'•READMN 


247 


7 




248 


STO 


00 


249 


»»NUM 


ROWS 


250 


XEQ 


'•IN 


251 


XEQ 


"NXT 


252 


XEQ 


"N XT 


253 


"NUM 


COLS 


254 


XEQ 


"IN 


255 


XEQ 


"NX? 


256 


RCL 


10 


257 


* 




258 


RCL 


04 


259 


+ 




260 


1 




261 


SU- 


00 


262 


RD N 




263 


STO 


IND 00 


264 


RCL 


09 


265 


+ 




266 


XEQ 


"SIZE? 


267 


FC?C 


25 


268 


PROMPT 


269 


RTN 
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HP 4 1 C SOURCE CODE 



LINEAR PROGRAMMING 



( 



SIZE? 



270 LBL "SIZE? 

271 "S I ZE> = 

272 ARCL X 

273 1 

274 - 

275 SF 25 

276 RCL IND X 

277 RTN 



"NXT 



278 


LBL 


"NXT 


279 


1 




280 


ST + 


00 


281 


+ 




282 


STO 


IND 00 


283 


RTN 





SETL 



284 


LBL "SETL 


285 


IE 03 


286 


/ 


287 


1 


288 


+ 


289 


STO IND Y 


290 


RTN 
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HP41C SOURCE CODE 



LINEAR PROGRAMMING 



READ JB 



291 


LBL 


"READJB 


292 


1 




293 


RCL 


08 


294 


XEQ 


"SETL 


295 


RCL 


13 


296 


STO 


00 


297 


FIX 


0 


298 


LBL 


01 


299 


"BASIC " 


300 


ARCL 


01 


301 


»!_, II 




302 


XEQ 


"IN 


303 


IS G 


01 


3 04 


GTO 


01 


305 


FIX 


4 


306 


RTN 
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HP41C SOURCE CODE: 



LINEAR PROGRAMMING 



( 

READA 



307 


LBL 


"READA 


308 


SF 10 


309 


1 




310 


RCL 


09 


311 


XEQ 


"SETL 


312 


LBL 


11 


3 13 


2 




314 


RCL 


12 


3 15 


XEQ 


"SETL 


316 


LBL 


12 


3 17 


FIX 


0 


318 


i» ip 




319 


ARCL 


01 


320 

321 


ARCL 


02 


322 


FIX 


4 


323 


LBL 


13 


324 


RCL 


01 


325 


1 




326 


— 




327 


RCL 


12 


328 


* 




329 


RCL 


02 


330 


+ 




331 


RCL 


04 


332 


+ 




333 


1 




334 


— 




335 


STO 


00 


3 36 


LBL 


14 


337 


FC? 


10 


338 


GTO 


16 


339 


XEQ 


"IN 


340 


GTO 


17 


341 


LBL 


16 


342 






343 


ARCL 


, IND 0 


344 


AVIEW 


345 


LBL 


17 


346 


IS G 


02 


347 


GTO 


12 


348 


IS G 


01 


349 


GTO 


1 1 


350 


CF 10 


351 


0 




352 


RTN 
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HP41C SOURCE CODE 



LINEAR PROGRAMMING 



( 



IN 



353 LBL "IN 

354 CF 22 

355 1 

356 ST+ 00 

357 RCL IND 00 

358 

359 PROMPT 

360 STO IND 00 

361 RTN 



INIT 



362 


LBL 


"INIT 


363 


CLRG 


364 


CF 


29 


365 


IE- 


04 


366 


STO 


07 


367 


19 




368 


STO 


04 


369 


RTN 





ERR 1 



370 


LBL 


"ERR1 


371 


0 




372 


LN 




373 


XEQ 


"READJB 



t 
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HP41C SOURCE CODS : 



LINEAR PROGRAMMING 







( 


CHECK 


374 


LE 


1l 


"CHECK 


375 


SF 


1 1 


376 


1 






377 


RCL 


12 


378 


XEO 


"SETL 


379 


RCL 


04 


380 


RCL 


09 


381 


RCL 


12 


382 


* 






383 








384 


STO 


14 


385 


STO 


15 


386 


LBL 


91 


3 87 


RCL 


14 


388 


RCL 


01 


389 


+ 






390 


0 






391 


STO 


IND Y 


392 


ISG 


01 


3 93 


GTO 


91 


394 


1 






3 95 


RCL 


08 


396 


XEQ 


"SETL 


3 97 


LBL 


92 


398 


CF 


07 


399 


RCL 


13 


400 


RCL 


01 


401 


♦ 






402 


RCL 


IND X 


403 


X<0? 


404 


SF 


07 


405 


ABS 




406 


STO 


00 


4 07 


X= 


= 0 


• 


408 


GTO 


"ERR 1 


409 


RCL 


12 


4 10 


X< = 


Y? 


411 


GTO 


"ERR 1 


412 


2 






413 


RCL 


09 


414 


XEQ 


"SETL 


415 


LBL 


93 


416 


CF 


08 


417 


RCL 


01 


418 


INT 




419 


RCL 


02 


420 


INT 




421 


X= 


: Y 


■7 


422 


SF 


08 


423 


1 






424 


— 






425 


RCL 


12 


4 26 


* 






427 


RCL 


00 
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HP 4 1 C SOURCE CODE 



LINEAR PROGRAMMING 



428 


+ 






429 


RCL 


04 




4 30 


+ 






431 


FC? 


08 




432 


0 






433 


FS? 


08 




434 


1 






435 


STO 


IND 


Y 


4 36 


IS G 


02 




437 


GTO 


93 




4 38 


FC? 


07 




439 


GTO 


98 




4 40 


CF 11 




441 


RCL 


14 




442 


RCL 


00 




443 


+ 






444 


1 






445 


STO 


IND 


Y 


446 


2 






4 47 


RCL 


12 




448 


XEQ 


"SETL 


449 


LBL 


96 




450 


RCL 


01 




451 


INT 






452 


1 






453 


- 






454 


RCL 


12 




455 


* 






456 


RCL 


02 




457 


+ 






458 


RCL 


04 




459 


+ 






460 


RCL 


IND 


X 


461 


CHS 






462 


RCL 


14 




463 


RCL 


02 




464 


+ 






465 


XOY 




466 


RCL 


IND 


Y 


467 


+ 






468 


STO 


IND 


Y 


469 


ISG 


02 




470 


GTO 


96 




471 


LBL 


98 




472 


ISG 


01 




473 


GTO 


92 




474 


FC? 


11 




475 


RTN 






476 


RCL 


04 




477 


RCL 


08 




478 


RCL 


12 




479 


* 






480 


+ 






481 


STO 


14 




482 


RTN 
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HP 4 1 C SOURCE CODE 



LINEAR PROGRAMMING 



i’ THE FOLLOWING TABLE DESCRIBES THE KEY REGISTER AND FLAG 

I ASSIGNMENTS MADE BY THIS PROGRAM: 

R00 = TEMPORARY REGISTER. HOLDS RECIPROCAL OF PIVOT 
( ELEMENT IN SUBROUTINE PIVOT. 

( R 0 1 = LOOP INDEX COUNTER 

( RO 2 = LOOP INDEX COUNTER 

{ R 0 3 = TEMPORARY REGISTER. HOLDS MIN VALUE IN FINDQ; 

? MAX VALUE IN FINDP; AND IS AN INTERMEDIATE 

1 ADDRESS CALCULATION VALUE IN PIVOT, 

r R04 = BASE REGISTER FOR INDIRECT ADDRESSING 

r (SET = 19) 

? RO 5 = THE PIVOT COLUMN NUMBER 

( R06 = THE PIVOT ROW NUMBER 

( R07 = EFFECTIVE ZERO LEVEL 

( R08 = M = NUMBER OF ROWS 

( RO 9 = M PLUS 1 

( RIO = M PLUS 2 

( R 1 1 = N = NUMBER OF VARIABLES 

R 1 2 = N PLUS 1 

( R 1 3 = BASE REGISTER FOR THE LOCATION OF THE 

( VECTOR JB WHICH CONTAINS POINTERS TO 

( WHICH VARIABLE IS BASIC IN EACH ROW. 

i 1 R 1 4 = ROW NUMBER OF THE OBJECTIVE FUNCTION; 

SET TO M PLUS 1 OR M PLUS 2 AS DETERMINED 
BY NEED FOR PHASE I. 

( R 1 5 = 3ASE REGISTER FOR THE LOCATION OF THE 

i‘ PHASE I OBJECTIVE FUNCTION, IF NEEDED. 

( R 1 6 = TEMPORARY REGISTER. 

R 1 7 = NUMBER OF PIVOTS PERFORMED. 

( R 1 8 = RESERVED FOR FUTURE USE. 

' R19-R177 = DATA STORAGE REGISTERS FOR ELEMENTS OF 

T THE TABLEAU AND THE JB VECTOR. 

i; FLAGS 

r F01 - F0 3 = SUBROUTINE EXECUTION FLAGS. BECAUSE THESE 
( FLAGS ARE VISIBLE IN THE DISPLAY THEY CAN 

( BE SET WHEN ENTERING A MAJOR SUBROUTINE AND 

r CLEARED WHEN LEAVING — GIVING THE USER A 

( VISUAL INDICATION OF WHAT IS HAPPENING 

INSIDE THE CALCULATOR. 

i- F01 — SUBROUTINE FINDQ 

( F02 — SUBROUTINE FINDP 

; F03— SUBROUTINE PIVOT 

’ F04 = WHEN SET, STOPS CALCULATOR AFTER EACH PIVOT. 

( F07 = USED AS TEMPORARY FLAG IN PIVOT AND CHECK ROUTINES. 

( F 1 0 = USED AS A TEMPORARY FLAG IN READ ROUTINES. 

( F 1 1 = WHEN SET, INDICATES PHASE II IS IN PROGRESS. 

( F29 = CONTROLS FORMAT OF DISPLAY SEPARATOR. 

483 END 
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APPENDIX C 

SUBROUTINES FOE READ ONLY MEMORY 

INTRODUCTION: 

The calculator subroutines described in this appendix 
perform functions which are frequently required by applica- 
tion programs and are therefore ideal candidates for use in 
a read only memory (ROM.) These routines were written espe- 
cially to illustrate the differences between read only mem- 
ory routines and routines designed for individual use via 
bar code or magnetic cards. These differences include more 
attention to entry and exit point options, an attempt to 
keep the size of the routines as compact as possible, and an 
attempt not to disturb the register stack if at all 
possible. 

These common subroutines are provided separately from 
application programs because when more than one application 
program uses the routines, as is recommended for these pro- 
grams, the use of a separate block of common functions saves 
space in the ROM overall. Also, by providing a convenient 
set of "macro” instructions, application programs can be 
constructed more quickly and easily. Because these 
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subroutines are used frequently, they have been individually 
optimized and tested to save memory space and execution 
time. 3y using these "macros" within an application pro- 
gram, the application programmer can be reasonably certain 
of their efficiency and reliability. 

Almost all user /calculator interface is handled by these 
routines. There is one set of subroutines which assumes the 
user has a printer, and one set which does not. Printer 
instructions are preceeded in the listings shown in the 
appendix by an (PRT: label. When not using these routines 
on read only memory, the user will load one set or the other 
(but not both), as appropriate for his/her application. In 
so doing, the user with the printer gets full benefit from 
it while the user without the printer pays no penalty in 
execution time or memory space for the calculator's print 
instructions. Also, to change from use of the printer to 
use of the calculator without it, the user needs only to 
read in the new common block — the application programs are 
retained in memory unchanged. The subroutines appear the 
same to any application program — giving the added benefit 
that any application program which uses them for input or 
output operations will automatically make good use of 
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the printer even if written by a programmer who did not 
explicitly consider a printer requirement. 

When using these common subroutines, a discipline is 
enforced upon the application program concerning use of the 
calculator memory registers. This saves the programmer 
from having to plan his "register map" from scratch for each 
new program. Also, it insures comparability across differ- 
ent application programs for similar data objects such as 
matrices and loop index counters. One of the most annoying 
problems with read only memory programs available from the 
calculator manufacturer is this lack of cross program com- 
patability. Conflict in the use of memory registers is the 
rule, rather than the exception; and it is frequently impos- 
sible to efficiently use more than one read only memory pro- 
gram as a subroutine within a user written program. A third 
reason why register assignment standards are advantageous is 
that they make it easier for the user of the calculator to 
remember the key register assignments and, if necessary, 
recall their contents during the execution cf an application 
program . 

Another function performed by this set cf common subrou- 
tines is to simplify the use of indirect addressing — a 
critical goal on the HPU1CV. 



109 



3ecause the common subroutines listed in this appendix 
are always called by application programs and never from the 
keyboard by the user, the typical user instructions are 
inappropriate. Instead, for the benefit of application 
programmers wishing to use the routines, the basic functions 
and structure of each are explained in subsequent sections 
of this appendix. 

Subroutine "IN” 

Subroutine "IN" is used as a general input and output 
interface between the user and the calculator. This subrou- 
tine has three alternative entry points which when called 
affect functions as follows: 

IN — Input mode (displays a question mark query) 

10 — Output mode (displays labeled data value) 

IX — Direct mode (input of value in x register) 

In particular, one entry point, "IN", may be called whenever 
an application program must query the user for a numeric 
input value. As such, it is a direct replacement for the 
PROMPT instruction organic to the calculator, but automati- 
cally prompts, verifies and stores the received value using 
an indirect address contained in register 00. The printer 
version of the subroutine will automatically label and print 
the final, verified data value recorded. 
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Subroutine "IN" uses register 00 as a data location 
pointer and automatically increments this register so that 
subseguent calls to the subroutine will result in sequential 
data manipulation. The application programmer must insure 
that register 00 contains a number equal to the storage reg- 
ister number prior to calling the subroutine. For example, 
if register 00 contains 17, "IN" will store the data in reg- 
ister 17. One of the major advantages of this routine is 
that the same subroutine may be used to verify and/or change 
the data previously recorded. Thus, separate edit routines 
are usually unnecessary. Pressing the R/S key without 
touching any other key leaves the va*lue stored unchanged. 
Pressing "1" and "+" and then "R/S" adds one to the current 
stored value, and so on for other function keys. Entering a 
new string of digits results in that new value being stored. 

An additional feature of this subroutine is the "verify" 
mode indicated by flag 05. Flag 05 is reserved for this 
purpose and is set "on" by a call to subroutine "VR"-- 
ancther of the subroutines in this common set. The verify 
mode is intended for use by a novice or other user who 
wishes to verify every data value as it is keyed into the 
calculator. The advantage is increased accuracy and 
confidence in the result. 
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Subroutine "D2" 

Subroutine "D2" is used to set up the index register for 
a program loop. This subroutine has two alternative entry 
points which when called, increment different index regis- 
ters as follows: 

D2 — Establishes register 02 as the index 
D 1 — Establishes register 01 as the index 

This subroutine is intended for use with the "ISG" loop 

structure which has the effect most like that of a FORTRAN 

"DO LOOP." For example, to execute a loop 20 times: 

• • • 

• • • 

20 

XEQ " D 1 
LBL 00 

■*« 

• • • • 

• • • • 

ISG 01 
GTO 00 
• • • 

• • • 

The advantage of this form of loop structure is that 
register 00 may be used within the loop for address calcula- 
tions. The first time the loop is addressed the integer 
portion of the number in register 00 will be 1, the second 
time it will be 2 and so on . There is no need to truncate 
the fractional portion of the number because the HP41C 
ignores the fractional component of a number when 
calculating a register address. Ose of index registers for 
address calcula-tion makes indirect addressing practical. 
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Registers 01 and 02 should be reserved for use as index 
registers by the application programmer. In most cases 
these two registers should prove sufficient. 



S ub routine " VR" 

Subroutine "VR” is used as a general purpose calculator 
initialization routine. This subroutine has three alterna- 
tive entry points which vary the amount of initialization 
performed as follows: 

VR — Sets the verify mode on, and the following: 

WR — Suppresses the audio tones, and the following: 

WS — Clears all memory. 

Sets the display for integers. 

Assigns statistical reaisers. 

Sets "zero" level for equality testing, and 
Sets base address for indirect addressing. 

In the printer version of this initialization routine, 
the subroutine prints a banner (usually the program name) 
which has been stored in the alpha register prior to calling 
the initialization subroutine. 

If flag 13 is set prior to calling the initialization 
routine, then the calling program must have placed the num- 
ber of data registers required to execute the program in the 
x register prior to calling the initialization routine. In 
this case, a check will automatically be performed using 
subroutine "SZ" described below. 

It is recommended that all application programs provide 
an alternate entry point which bypasses the initialization 
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szep. Then if a data error is encountered, the user may 
review the data entered into the calculator by simply press- 
ing the return key once after every prompt. This procedure 
works because subroutine "IN" recalls the stored value prior 
to prompting the user. When the user presses the clear key, 
the alphabetic prompt is removed and the existing data value 
revealed. 

Sub r out ine " SZ " 

Subroutine "SZ" is used to test if sufficient numbers of 
data registers are available to run an application program. 
This subroutine may be either called directly or as part of 
the initialization routine described above. 

Subroutine "SR" 

Subroutine "ER" is called whenever the application pro- 
gram encounters an error--usually in the input. data. A 
prompt is displayed and an audio tone sounded, provided flag 
26 (the audio enable flag) has not been cleared by the 
initialization routine described in paragraph D. 

Su brou t ine "SORT" 

This subroutine is included to illustrate the use of a 
stack register table in the program comments, but it also 
represents a useful utility routine. The sorting algorithm 
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used is the shell sort [Ref. 6: pp. 84-95] which gives 

reasonably fast sorting times with a very small program 
size. All conventions such as base register in R04 and num- 
ber of data points in R15 are assumed by this subroutine. A 
complete list of all such register assignments is listed at 
the end of the program listing. 

Subrout ines "PUT" and "GET" 

These two small routines provide a useful capability to 
store and recall up to three integers between 0 and 999 in 
one data register. This means that if you are manipulating 
a spread sheet of small, positive integers you can store the 
same data in one third the space. Of course, run time is 
degraded (about 20 seconds for every 100 data references.) 

To store a value, assuming the base register has been 
defined, just press: 

value ENTER | point-number XEQ "PUT 
To recall a value, simply press: 

point-number XEQ "GET 
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FP41C SOURCE CODE: 



COMMON SUBROUTINES 



( ( 
( 

( 

( 

( 

l 



j 



IN 





( 




— ) 


L 


LdL " IN 




(INPUT MODE — 01 SPLAY LABEL AND ? 


2 


SF 1C 




(SET QUERY ONCE FLAG 


3 


GTQ C5 






4 


LdL "IC 




(OUTPUT MODE — DISPLAY LABEL AND DATA 


5 


CF 10 




(INSURE NO QUERY 


6 


LdL 05 






7 


RCL I NO 00 




(ASSUMES ROD POINTS TO STORAGE REG 


e 


L3L "IX 




(DIRECT MODE — ASSUMES X REG FOLDS DATA 


9 


A3TC 05 




(ASSUMES LABEL SET UP BY CALLING PROG 


10 


L3L 01 






i i 


"-1 = 






12 


FS? 10 




(QUERY OR DISPLAY VALUE? 


13 


"-»? 






14 


FC? 10 






15 


ARCL INC 00 




(DISPLAY DATA VALUE 


16 


CF 22 




(DIGIT ENTRY DETECTION FLAG 




( P R T : FS? 


10 






1 PPT : CF 21 


(DISABLE PRINTER FOR QUERY 


17 


PROMPT 








I PRT : SF 21 




13 


ST C IND 00 




(STORE INPUT VALUE 


19 


CLA 




(PREPARE ALPHA REG FOR NEXT PROMPT 


20 


FC? 05 




(NOT VERIFY MODE? 


21 


GTC 03 








( FOLLOWING 


I S 


ACTION TAKEN WHEN IN VERIFY MODE 


22 


FS? 22 




(ANY INPUT DETECTED? 


23 


GTO 02 




(REDISPLAY VALUE IF USER INPUT 


24 


FC? 10 




(DISPLAY VALUE ONCE IF NO INPUT 


25 


GTC 04 






26 


LdL 02 




(DISPLAY LABEL AND DATA VALUE 


27 


CF 10 




(NOTE MUST HAVE QUERIED ONCE 


28 


ARC L 05 






29 


GTO 01 








( FCLLCW ING 


IS 


ACTION TAKEN FOR NON-VERIFY MODE 


30 


LB L 03 






2 1 


FS?C 10 




(INPUT MODE? 


22 


GTG 04 




(IF INPUT MODE, EXIT ROUTINE 


33 


FS? 22 




(WAS THERE INPUT ? 


34 


GTC 02 




(USER CHANGE C VALUE, SC VERIFY. 


35 


LdL 04 




(PREPARE TO EXIT 




( PR T : FS? 


55 


(PRINTER ATTACHED? 




( PRT : GTO 


0 6 






(FFT: LBL 


05 




36 


ISG 00 




(INCREMENT POINTER FOR NEXT IC OPERATION 


37 


RTN 






38 


RTN 




(FINAL VALUE IS IN X REG UPON EXIT 




(PRT: LBL 


06 


( PRINT THE FINAL VALUE 




(PRT: CLA 








(PRT: ARCL 


05 






(PRT: 








(PRT: ARCL 


IND 


00 




(PRT: PRA 








(PRT: GTO 


05 
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HP 4 1C 5CLRCE CODE: 



C CM MON SUBROUTINES 



( ( 



D 2 



39 LBL "02 

40 2 

41 GTC 07 

42 L3L "01 

43 1 

44 L3L 07 

45 XOY 

46 1E3 

47 / 

46 1 
49 + 

5 C STG INC Y 
51 RTN 



) 

) 

) 

) 

) 

) 

(SETUP LOOP USING REG 2 AS INDEX COUNTER 



(SETUP LOOP USING REG 1 AS INDEX COUNTER 



(NUMBER LOUP ITERATIONS MUST G5 IN X 
(PRIOR TO CALLING THIS SUBROUTINE. 
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HP 4 L C SOURCE CODE: 



COMMON SUBROUTINES 



{ ( 

( 

( 

( VR 

( 

{ 

( 

5 2 LB L "VR 

53 SF C 5 

54 L3L "WR 



5 5 


CF 26 




56 


LBL "«S 




57 


CLRG 




58 


£R EG 10 




53 


CF 29 




60 


1E-C4 




61 


ST C 03 




62 


15 




63 


STO 04 




6 4 


1 




to 


+ 




66 


STO GO 






( PRT: 


AD V 




( PRT: 


SF 12 




( PRT: 


FC? 55 




( FRT: 


PRA 




(PRT : 


CF 12 




(PRT: 


ADV 


67 


FC? 13 




6 8 


RTN 




69 


X<> Z 





j 

) 

) 

) 

) 

j 

(SET VERIFY MODE ON 

(INITIALIZE FOR CLASSROOM USE — NO AUCIO 
(STANDARD INITIALIZATION ROUTINE FOLLOWS 



(SETS DISPLAY — NO DECIMAL POINT IF INTEGER 
(EFFECTIVE ZERO LEVEL 

(NORMAL INDIRECT ADDRESS BASE REGISTER. 

(this is i less than the number of the 



(FIRST REGISTER WHERE DATA IS STORED. 

( D0U8LE W IDE PRINTING 
(PRINTER ATTACHED? 

(SET BACK TO SINGLE WIDE PRINT 

(HAS INITIALIZATION ROUTINE BEEN ASKED TO 
(VERIFY SIZE? 
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FP4 LC GCU RC E COCE: 



COMMON SUBROUTINES 



( ( ) 

( ) 

( ) 

( sz 

( J 

( ) 

70 131 "SZ 

71 "SET SZ= 

72 ARCl X 

73 1 

74 - 

75 SF 25 (PREPARE TO IGNORE ERROR 

76 RCL IND X (TEST FLAG 25 TO SEE IF SUFFICIENT 

77 FC?C 25 

78 PRCPPT 

75 RTN (SIZE EXISTS. 

( j 

( ) 

( ER 

( ) 

( ) 

( j 

30 LdL "ER (DISPLAY "DATA ERROR" PROMPT C SOUND TONE. 

3 L 0 

c 2 LN (BEST *AY TO DETERMINE *iHERE ERROR 

S3 TONE 2 (OCCURED IS TO HIT THE SST KEY ONCE, 

34 RTN (THEN GO INTO PR GM MODE. 
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HP 4 L C SOURCE CGDE 



C CM MON SUBROUTINES 



l ) 

( ) 

( SORT 



35 


LB L 


"SCRT 
















6 b 


RCL 


15 






(RECALL 


NUMBER OF DATA POINTS 




67 


ST C 


01 






( DEF INE 


A = 


"MIDPOINT 


" CF NU- V 3ER POINTS 


ea 


L3L 


00 


















as 


RCL 


Cl 






(RECALL 


MIDPOINT 








so 


2 




















51 


/ 




















52 


INT 




















5 3 


STC 


01 






(A = INTIA/2 ) 








54 


X=C? 






( TEST TO SEE IF 


LIST 


SORTED 




C e, 


RTN 




















56 


1 




















57 


STG 


02 






(8=1 


— 


RESET 


LEFT 


SHELL BOUNDARY 


58 


LBL 


01 




( 


STACK 


TABLE FULLGwS: 






55 


STC 


03 




( 


C=B 




8 








ICC 


LBL 


02 




( 


C 




B 








1CI 


RCL 


01 




( 


A 




C 




3 




1C2 


+ 






{ 


D = C+A 




B 








IC3 


RCL 


04 




( 


BASE 




D 




B 




104 


+ 






( 


A COR 0 




B 








ICS 


RCL 


IN0 


X 


( 


X (D ) 




A DDR 


D 


8 




IC6 


RCL 


03 




( 


C 




X (D) 




A DDR D 3 




1C 7 


RCL 


04 




( 


BASE 




C 




X ( D ) ADDR 


D 


ica 


+ 






( 


A DDR C 




X (D) 




A DDR D 




IC5 


XOY 




( 


X (D) 




A DDR 


r 


ADDR D 




l 10 


RCL 


INC 


Y 


( 


X (C ) 




X ( D) 




ADDR C ADDR 


D 


111 


X< = 


Y? 


















112 


GTQ 


03 




( FOLLOWING 


WILL 


. INTERCHANGE X(C) AND X ( D ) 


1 13 


STC 


I NO 


T 


( 


X (C ) 




X (D) 




ADDR C ADDR 


D 
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XOY 




( 


X (D) 




X (CJ 




ADDR C ADDR 


D 


1 15 


STG 


INC 


l 


( 


X (0 ) 




X (C) 




ADDR C ADDR 


D 
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RCL 


0 3 




( 


C 












1 17 


RCL 


Cl 




( 


A 




C 
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— 






( 


C =C - A 












115 


STC 


03 




( 


C 












1 2 C 


X>C 


t 


















121 


GTC 


02 


















122 


LBL 


03 


















123 


RCL 


15 




( 


N 












124 


RCL 


01 




( 


A 




N 








125 


— 






< 


E=N -A 












126 


RCL 


02 




( 


B 




E 








127 


1 






( 


1 




B 




E 




128 


+ 






( 


3+1 




E 








125 


STC 


02 




( 


ri =3 + 1 




E 








130 


X< = 


Y? 


















131 


GTC 


01 


















132 


GTG 


CC 




















( 






SCRT 


USES 


THE 


FCLLCWI NG 


REGISTERS: 






( 








R 01 — 


A 












( 








RC2 — 


B 












( 








R 03 — 


C 












( 








R 04 — 


IND I RECT 


ADDRESS BASE 






( 








R 15 — 


NUMBER DATA POINTS (SET 3Y 


£+ ) 
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EP41C 3CIRCE CQDt: 



COMMON SUBROUTINES 



/ 



< ( 

( 

( 

( SC 

( 



( 

{ 



123 


LB t 


"SC 


124 


RCL 


00 


125 


1 




126 


+ 




127 


ST C 


05 


128 


-1EC3 


139 


STC 


10 


140 


-1 




141 


3 




142 


/ 




143 


STG 


09 



) 

) 

) 

) 

) 

) 



121 



b F 4 1C SCI PC E CODE: CCMMUN 

( ( 

( J 

( 

( per 

( 

( 



144 


LB L 


"PUT 


145 


"VALUE 


146 

147 


11 

STC 


00 


146 


XEQ 


"I A 


149 


SF 


C2 




GET 



( 

( 

( 

150 LB L "GET 

151 "BIT REG 

152 PRC V P T 



SU BR CUT I NE3 
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HP 4 1 C SCLBCE CODE: 



COMMON SUBROUTINES 



( ( 






SA 






15 3 


LB L 


"SA 


154 


■7 




155 


7 




156 


INT 




157 


STC 


07 


158 


LASIX 


l 59 


— 




16C 


RCL 


09 


161 


/ 




162 


12 




163 


+ 




164 


STO 


C 3 


165 


RCL 


05 


16 6 


ST + 


07 


167 


RCL 


INQ 


168 


INT 




169 


STC 


14 


1 1C 


LAS 


IX 


171 


— 




1 72 


RCL 


10 


173 


* 




174 


INT 




175 


STO 


13 


176 


LAS 


IX 


177 


— 




178 


RCL 


10 


179 


* 




iao 


INT 




181 


STC 


12 


182 


RCL 


11 


183 


RCL 


INC 


1 84 


STO 


11 


165 


FC?C 02 


186 


RTN 




187 


XOT 


188 


STC 


IND 


189 


RCL 


12 


ISC 


RCL 


10 


191 


CHS 




192 


/ 




1 S3 


RCL 


13 


IS4 


+ 




195 


RCL 


10 


1 S 6 


CHS 




197 


/ 




198 


RCL 


14 


l S 9 


♦ 




2 CO 


STC 


INC 


2 C 1 


RTN 
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HP41C SOURCE CODE: 



COMMON SUBROUTINES 



202 



( ( 

( 

( 

I 
( 

I 
( 

{ 

( 
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( 

( 

( 
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( 

( 

t 

( 

( 

( 

( 

( 

( 

( 

{ 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

( 

END 



THE FOLLOWING TABLE 
ASSIGNMENTS MADE BY 



DESCRIBES THE 
THIS PROGRAM: 



KEY REGISTER AND FLAG 



RGO = INDIRECT ADDRESS FOR STGRAGE OF INPUT DATA 

R 01 = LOCP index counter 

R02 = LCCP INDEX COUNTER 

R03 = EFFECTIVE ZERO LEVEL — USE AS TEMP IF NA 
R04 = dASE REGISTER FOR INDIRECT ADDRESSING (15) 

R05 = TEMP REGISTER FOR ALPHA PROMPT 
R 06 - R09 = APPLICATION PROGRAM TEMP REGISTERS 
RIO - A 1 5 = STATISTICAL REG I STERS--U SE AS TEMP IF NA 
R 16 . . . . = STORAGE OF DATA VIA INDIRECT ADDRESSING 



FCO-FOA = 



FLAGS 

SUBROUTINE EXECUTION 
FLAGS ARE VI SI BLE IN 
BE SET WHEN ENTERING 
CLEARED wHEN LEAVING 
VISUAL INDICATION OF 
INSI DE THE CALCULATOR 
FLAGS IF THIS IS NOT REQUIRED. 



FLAGS. BECAUSE THESE 
THE DISPLAY THEY CAN 
A MAJOR SUBROUTINE AND 
— GIVING T (E USER A 
WHAT IS HAPPENING 
USE AS TEMPORARY 



F05 = VERIFY INPUT MODE. ''ON'* WHEN SET. WHEN SET 

AFTER EVERY DATA VALUE IS ENTERED, THE CALC. 
WILL ECHO THE PROMPT AND DATA VALUE ENTERED. 
IF VALUE IS CORRECT, SIMPLY PRESS R/S KEY, 
OTHERWISE ENTER A CORRECTED VALUE AND 
CALCULATOR WILL AGAIN ASK FOR VER I F I CAT I ON . 



F 10 = 

FI 1 = 
F 12 = 
F 1 3 = 



F 14 - 



1C 



INDICATES 



USED AS A TEMPORARY FLAG INSIDE 
10 QUERY PROMPT IS DESIRED. 

AUTOMATIC EXECUTION FLAG -- DON’T USE EVER 
J0U3LE WIDE PRINTING -- LOCAL WITHIN "IO”. 

WHEN SET MEANS MAIN ROUTINE ASKING "VR" FOR 
AUTOMATIC SIZE CHECK AFTER INITIALIZATION. 
APPLICATION PROGRAM MAY USE FREELY AFTER INIT. 
F20 = AVAILABLE FOR APPLICATION PROGRAM USE. 
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APPENDIX D 

THE CROSS COMPILER PROGRAM AND COMMAND PROCESSOR 
DESIGN METHODOLOGY: 

This appendix discusses the design methodology used dur- 
ing construction of both the cross compiler program and the 
command processor, which is an IBM EXEC II program which 
provides an interactive programming environment for users of 
the system. 

Blazie's compiler for the HP65 calculator [Ref. 9] rep- 
resents one of the first attempts to provide a compiler for 
calculator programs. Both Carvalho [Ref. 10: pp. 25-29] 

and McNeal [Ref. 11: pp. 148-178] have published BASIC lan- 

guage programs which cross compile HP41C7 ins-ructions on a 
microcomputer for output to a line printer which can print 
acceptable bar code. While these referenced programs pro- 
vided valuable insights into the problem, especially into 
the special character istics of the HP41CV instruction set, 
none was exactly suited to the needs of this study. Because 
the Versatec plotter at the Naval Postgraduate School could 
be easily used only by FORTRAN programs, FORTRAN seemed the 
computer language of choice for this project. Both programs 
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were written with limited objectives and neither would have 
easily supported the extensions desired. Extensions planned 
for implementation included: 

• An extended instruction set. 

• In code comments. 

• Extensive error checking. 

• Compatability with the Emulator. 

• Synthetic Instructions [Ref. 1]. 

• Instruction macro's. 

Having decided to code an original cross compiler, a 
design methodology which would capitalize on the advantages 
of FORTRAN was planned. FORTRAN'S major deficiency for use 
in constructing a compiler of any type is its lack of 
alpha-numeric string handling capabilities. Rather than 
struggle with the lack of string functions, it was decided 
to code the necessary string functions as separate subrou- 
tines. This decision was reinforced countless times 
throughout the process of writing the compiler. The string 
function subroutines have been used not only in the cross 
compiler, but in many other FORTRAN programs since they were 
originally written. In fact, many persons who have no 



129 



interest in the HP41CV cross compiler may find the set of 
string functions listed in this thesis to be a valuable set 
of utility routines to be used to augment FORTRAM. The gen- 
eral convention used througout the string function subrou- 
tines is that an alphabetic string may be represented as a 
vector of two byte integer variables used to store the char- 
acters and a single four byte integer variable used to 
store the length of the string. 

One of the major advantages of the cross compiler is its 
ability to handle comments integrated within the HP41CV 
source code. This feature is critical to the clarification 
of the logical structure of the HP41CV programs. Because 
the parenthesis is nor a valid HP41CV character, it was cho- 
sen as the comment indicator character. A comment may occur 
beginning at the first column on an input line or anywhere 
after an HP41CV instruction. The comment must follow the 
instruction because everything after the comment mark out to 
the end of the input line is considered part of the comment. 

The control the user has over the output listing is also 
one of the advantages of the cross compiler. When twc com- 
ment indicator marks are placed in positions one and two of 
the input line, the compiler will force a page eject when 
printing the output listing. In addition, the user can vary 
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the number of output lines per page and cause useful tanners 
to be placed adjacent to program labels for ease of 
recognitio n. 

Altogether there are twenty-four subroutines and a main 
program which constitute the cross compiler. The source 
code for each of these routines is provided in the second 
section of this appendix. Each subroutine begins with a 
statement concerning its function and construction. Accord- 
ingly, no general description of each subroutine will be 
repeated here. However, subroutine COMP deserves special 
attention, for it is the master lexicographic analyzer for 
the compiler and would also interface the user with the emu- 
lator. Its function is to receive a single line of HP41CV 
source code and identify it. COMP considers all HP41CV 
instructions to be of one of three types. The first cat- 
egory are the single byte instructions with no operands that 
can be compiled by a simple table look up. COMP has been 
constructed so that the instruction set can be extended at 
any time simply by increasing the size of this table. In 
this way abbreviated or altered command names could be eas- 
ily used. The second category of instructions are the 
multi-byte instructions which require a table lookup and the 
translation of one or more operands, including possibly an 



131 



indirect instruction indicator. Ihe table examined by the 
compiler is the same as for the category one instructions, 
and a code is given in the table which indicates to the 
compiler the number of operands which are required with each 
instruction. A syntax check is then made in subroutines 
10 NE and ITWO to insure that the number and characteris- 
tics of each operand are appropriate for the given instruc- 
tion. One of the major advantages of the use of the cross 
compiler is the syntax and error checking that is performed 
during the compilation process. The third type of instruc- 
tion represents the exceptional instructions that are so 
difficult to compile that they require separate subroutines 
for efficient compilation. These instructions include stor- 
age and recall of data, program labels and program flew con- 
trol statements such as goto and execute. 

In order to provide an efficient programming command 
system for the compiler that would minimize the need to know 
technical details about the operation of the compiler, an 
IBM EXEC II program was written. This program not only 
interfaces the user to the compiler, bur it also provides on 
line user instructions as t o how to use the system. 

Included in this command processor is a command menu which 
gives the format and short description for each command. 
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Another command, help, provides more detailed information 
about each command. When a novice user first enters the 
exec, or types the name of the exec followed by a question 
mark, then he receives a four page narative description of 
what the system is, how it works, and what actions he must 
take to write a successful HP41CV application program. 
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THE ROLL C * ING IS THE DATA SET RE AD BY THE CROSS COMPILER 
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64 


- 


65 


* 


66 


/ 
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X< Y 7 


68 


X > Y ? 


69 
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LN 
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29 3 1 

99 4 1 

97 4 1 

75 4 1 

91 5 1 

5 7 4 1 

74 4 1 

20 5 1 
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33 3 1 

63 3 1 

46 4 1 
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25 3 1 

16 4 1 

64 3 1 
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14 3 1 

60 3 1 
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12 3 1 

69 3 2 

33 3 2 

79 3 2 
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80 3 2 

82 3 2 
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34 3 2 
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13 4 2 

15 4 2 

42 3 

73 3 

36 3 

86 4 
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